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EXECUTIVE  SUMMARY 


This  document  presents  a  final  report  to  the  Office  of  Naval  Research  for  the  research  program  enti¬ 
tled  "A  Theoretical  Search  For  Supen^elocity  Semiconductors".  This  program  has  been  funded  by  ONR 
since  1974  in  the  Department  of  Electrical  and  Computer  Engineering  at  N.C.  State  UniversityTThe 
research  has  resulted  in  more  than  65  refereed  publications  and  numerous  conference  presentations  from 


its  inceptions  Major  contributions  to  the  field  of  hot  electron  transport  and  semiconductor  device  model¬ 
ing  have  been  achieved,  new  computational  methods  have  been  developed  (e.g.  path  integral  Monte  Carlo 
techniques),  and  the  work  has  helped  stimulate  commercial  ventures  in  the  applications  of  quaternary 
semiconductor  materials  to  electronic  and  optical  devices;)  In  addition,  there  have  been  twenty-four  Ph.D. 


and  M.S.  students  who  have  received  degrees  at  N;C.  State  University  with  research  support  from  this 


contract.  Three  visiting  faculty  members.-from  Japan  came  to  the  University  to  work  with  the  faculty 
investigators  supported  imder  this-ONR  contract  during  the  1979-1983  time  period.  A  visiting  professor 
from  the  French  CNRS  Microstructures  and  Microelectronics  Laboratory  in  Bagneux  (near  Paris)  spent  a 
sabbatical  year  at- N.C.  State  during  1988-89,  and  he  devoted  full-time  working  on  this  program  at  no  cost 
to  ONR.  ^During  the  current  funding  period,  a  visiting  scholar  from  China  is  a  member  of  our  research 

group  working  on  projects  which  directly  impact  this  ONR  program. 

/ 

This  initial  phases  of  this  work  centered  around  the  development  of  Monte  Carlo  simulation  tech¬ 
niques  which  allow  the  study  of  detailed  physics  of  hot  electron  transport  in  a  variety  of  compound  sem¬ 
iconductor  materials,  -^^e  original  emphases  were  concerned  with  electronic  materials  phenomena.  Later 
work  considered  the  utilization  of  these  materials  in  realistic  device  structures  where  physical  boundary 
conditions  must  be  imposed  on  the  carrier  transport.  More  recently,  the  work  has  focused  on  the  domain 
of  ultra-small  materials  and  device  phenomena  where  microscopic  non-local  transient  effects  such  as 
velocity  overshoot,  ballistic  and  nearly-ballistic  transport,  and  quantum  transport  become  important  or 
dominant.  During  the  past  four  years  we  have  begun  new  research  into  the  applications  of  the  Feynman 
"integral  over  paths"  approach  to  quantum  transport  as  well  as  the  study  of  hot  electron  effects  in  new 


device  structures,  such  as  the  hot  electron  spectrometer,  heterojunction  bipolar  transistor,  small  dimension 
metal-semiconductor-metal  photodetectors,  and  delta-doped  high  electron  mobility  transistors.  Also,  we 
have  been  exploring  some  new  approaches  to  device  modeling  which  combine  the  Monte  Carlo  method 
with  the  method  of  moments  of  the  Boltzmann  transport  equation  (hydrodynamic  transport  model)  for 
studying  specific  device  structmes,  such  as  small-dimension  n'*'  -n-n'^  majority  carrier  devices  and  the  high 
electron  mobility  transistor.  During  the  past  year  we  have  incorporated  quantum  correction  terms  into  the 
hydrodynamic  model  and  applied  this  model  to  resonant  tunneling  structures.  Quite  recently,  we  have 
applied  a  new  ansatz  distribution  function  as  a  constimtive  relation  to  close  the  moment  equations  in  the 
hydrodynamic  transport  model.  Initial  results  of  this  approach  have  been  physically  satisfying  and  com¬ 
putationally  promising.  New  work  on  the  theory  of  optical  phonon  modes  in  heterostructures  has  demon¬ 
strated  the  importance  of  the  effects  of  reduced  dimensionality  on  transport  dominated  by  the  LO  phonon 
interaction  near  heterointerfaces.  As  a  result  of  this  work,  a  new  concept  for  velocity  enhancement  i# 
pseudomorphic  device  structmes  has  been  presented  to  ONR. 

During  the  last  report  period,^we  focused  our  efforts  on  new  research  in  four  areas  along  with  the 
continuation  of  our  basic  research  directions  encompassing  the  study  of  hot  electron  transport  in  materials 
and  devices.  These  four  areas  included  I)  the  study  of  quantum  transport  in  resonant  tunneling  and 
related  ultra-small  heterobarrier  device  structures  using  the  Feynman  path  integral  and  the  transmission 
matrix  method,  2)  the  study  of  the  physics  of  hot  electron  and  heterojunction  device  structures,  3)  some 
new  approaches  to  the  merging  of  Monte  Carlo  methods  with  moment  equation  methods  with  resulting 


improvements  to  the  hydrodynamic  transport  model,  and  4)  new  Monte  Carlo  simulation  in  order  to  study 
pseudomorphic  devices,  delta-doped  structures,  and  submicron  MOSFETs.  )  c? - 
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1.0  INTRODUCTION 

In  October,  1974,  the  Office  of  Naval  Research  initiated  sponsorship  of  a  basic  research  program  in 
the  Department  of  Electrical  and  Computer  Engmeering  at  North  Carolina  State  University.  The  general 
goal  of  this  research  program  has  been  and  continues  to  be  the  investigation  of  high-speed  carrier  tran¬ 
sport  in  III-V  compound  semiconductors,  Ill-V  alloy  materials,  other  advanced  electronic  and  optical 
materials,  and  novel  device  slructmes  which  utilize  these  materials.  Four  faculty  members  at  N.  C.  State 
have  been  primarily  involved  and  supported  by  this  project.  They  are  M.  A.  Littlejohn,  J.  R.  IRiuser  and 
T.  H.  Glisson,  presently  Professors  in  the  ECE  Department,  and  Dr.  K.  W.  Kim  who  joined  the  ECE 
faculty  in  August,  1988  and  who  has  served  as  a  PI  on  the  effort  since  January,  1989.  This  research  pro¬ 
gram  has  made  significant  contributions  to  the  imderstanding  and  knowledge  base  of  hot  electron  tran- 

• 

sport  in  materials  and  devices.  It  has  provided  scientific  guidance  to  the  U.  S.  Navy  in  the  formulation  of 
a  part  of  its  basic  and  applied  research  program.  Numerical  concepts  developed  under  this  project  have 
been  transferred  to  other  Universities,  including  the  University  of  Illinois,  and  we  continue  these  inter- 
institutional  collaborations.  During  the  past  year,  two  faculty  members  and  two  graduate  students  from 
NCSU  spent  extended  time  with  Prof.  K.  Hess.  Director  of  the  National  Center  for  Computational  Elec¬ 
tronics  at  the  University  of  Illinois.  Also,  during  the  last  year,  collaborations  with  Dr.  H.  L.  Grubin  of 
Scientific  Research  As.sociates.  Inc.  and  Dr.  G.  J.  lafrate  of  the  U.  S.  Army  Research  Office  in  the  area  of 
quantum  transport  in  semiconductor  devices  have  been  initiated.  Also,  our  research  results  have  helped 
stimulate  commercial  venhues,  particularly  in  the  development  of  GalnAsP-based  materials  and  devices. 
To-date,  this  program  has  resulted  in  57  referred  publications  in  the  literature,  3  additional  manuscripts 
are  currently  in  press,  and  8  manuscripts  have  been  submitted  or  are  in  preparation  for  submission  to  the 
technical  literature.  A  listing  of  these  publications  is  given  in  Appendix  A.  In  addition,  numerous  invited 
talks  and  presentations  have  been  given  at  conference'  and  workshops  throughout  the  United  States  and 
in  other  countries.  The  program  has  contributed  significantly  to  the  educational  program  at  N.  C.  State 
University  with  more  than  twenty-five  Ph.D.,  M.S.  and  undergraduate  students  having  received  support 
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under  this  ONR  contract.  Cuirently,  four  Ph.D.  students  are  working  toward  their  degrees  on  this  project. 
Two  of  these  Ph.D.  students  are  U.S.  Citizens.  A  visiting  scholar  from  China  is  also  working  on  research 
related  to  this  program. 

This  ONR  program  has  been  efficient  and  productive.  The  high  quality  of  this  research  will  con¬ 
tinue  to  be  maintained  in  the  future.  This  report  will  summarize  the  progress  and  accomplishments  made 
during  the  past  twelve-month  ONR  reporting  period. 

2.0  RESEARCH  RESULTS 

2.1  Background 

Since  its  inception,  semiconductor  technology  has  been  stimulated  by  requirements  for  electronic 

» 

systems  with  ever-increasing  capabilities  to  process  information  faster,  more  functionally  and  more 
efficiently.  These  requirements  have  motivated  the  scaling  down  of  integrated  circuit  (IC)  device  dimen¬ 
sions  into  the  submicron  (less  than  ten  thousand  angstroms)  and  ultrasubmicron  (less  than  one  thousand 
angstroms)  regions.  Today,  we  have  entered  an  era  where  nanostructure  physics  and  fabrication  motivate 
our  research  efforts  in  semiconductor  electronics  [1],  As  fabrication  technology  has  allowed  such  devices 
to  be  realized,  many  new  and  fundamental  questions  have  emerged  concerning  the  underlying  physics  of 
small  (atomic  level)  dimensions  in  semiconductor  devices.  Important  issues  now  under  consideration  for 
ultrasubmicron  devices  include  nonequilibrium  transport  dealing  with  such  topics  as  quasi-ballistic  tran¬ 
sport,  overshoot  phenomena  and  quantum  transport.  A  great  deal  of  progress  has  been  achieved  in  our 
understanding  of  these  important  device  effects,  although  major  work  remains  to  be  done  as  our  ability  to 
fabricate  very  small  electronic  device  structures  continues  to  expand  and  mature  [2,3]. 

The  ability  to  fabricate  small  devices  has  been  continually  refined  over  the  last  ten  years  through 
impressive  improvements  in  materials  growth  technologies.  Molecular  beam  epitaxy  (MBE),  metalor- 
ganic  chemical  vapor  deposition  (OMCVD),  and  atomic  layer  epitaxy  (ALE)  have  provided  the  capability 
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to  fabricate  a  wide  variety  of  materials  and  heterostructure  combinations  with  near  perfect  interfaces, 
doping  control  and  compositional  uniformity  witii  atomic  level  dimensions.  The  development  of  ALE 
may  very  well  prove  to  be  the  ultimate  growth  technology  since  it  allows  the  deposition  of  one  monolayer 
of  device  quality  materials  through  a  controllable,  self-limiting  mechanism,  and  is  especially  useful  for 
the  deposition  of  heterojunctions  [4j.  The  ability  to  grow  layers  with  dimensions  of  a  few  angstroms 
opens  the  domain  of  quantum  transport  to  experimental  study  and  verification.  Thus,  topics  resulting 
from  size  quantization  in  condensed  matter  must  be  investigated  from  theoretical  viewpoints  with  tools 
which  are  either  partially  developed  or  through  the  development  of  new  tools  which  are  not  now  avail¬ 
able.  Quantization  effects  arising  from  geometrical  size  constraints,  proximity  effects  resulting  from 
closely  packed  arrays  of  devices,  and  general  solid-state  considerations  not  heretofore  considered  ques¬ 
tionable  (effective  mass  approximation,  the  role  of  contacts  and  the  like)  must  be  addressed  from  a  funda¬ 
mental  point  of  view.  Moreover,  from  a  device  physics  point-of-view,  it  is  desirable  to  have  a  micros* 
copic  description  of  the  physics  of  small  dimensions  which  is  amenable  to  phenomenological  treatment, 
so  that  its  properties  can  be  meaningfully  incorporated  into  futuristic  device  concepts  and  simulations. 

Theoretical  methods  to  address  carrier  transport  have  also  progressed  rapidly  over  the  last  ten- 
fifteen  years,  in  a  manner  similar  to  research  in  semiconductor  thin  film  epitaxial  growth  technology.  In 
fact,  this  is  a  natural  progression  in  many  ways  and  is  to  be  expected.  The  progress  achieved  in  materials 
growth  of  structures  with  quantum  dimensions  dictates  that  new  approaches  be  developed  and  refined  to 
study  quantum  transport  phenomena.  However,  a  significant  change  in  direction  is  now  warranted.  Past 
transport  theory  and  device  modeling  approaches  have  relied  on  particle  or  quasi-particle  approaches 
where  the  electrons  are  treated  as  rigid  mobile  entities  which  undergo  interactions  with  the  transport 
medium.  Tlic  treatment  of  the  interaction  often  involves  wave  concepts.  However,  the  model  is  basically 
a  particle  model.  In  the  current  regime  of  quantum  transport,  we  may  no  longer  be  able  to  consider  the 
carriers  as  particles.  It  is  quite  likely  that  their  physical  behavior  will  be  governed  either  partially  or  com¬ 
pletely  by  wave  phenomena. 
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The  quasi-particle  methods  attempt  to  retain  as  much  of  the  classical  formalism  as  possible  in  order 
to  be  able  to  express  results  in  terms  of  parameters  which  are  of  the  greatest  experimental  interest,  such  as 
carrier  velocity  and  diffusion  constant.  This  is  a  flexible  approach.  However,  much  care  is  required  to 
ensure  that  all  important  effects  are  properly  included  because  of  the  approximations  involved  in  the  for¬ 
mulations.  On  the  other  hand,  more  fundamental  quantum  approaches,  such  as  operator-eigenfunction 
methods,  adhere  closely  to  the  actual  quantum  states  present  in  the  device  structure  when  scattering  is  not 
included.  Scattering  processes  (dissipation)  can  be  added  using  pertinbation  theory  from  quantum 
mechanics.  These  techniques  can  obtain  the  greatest  sensitivity  to  the  resulting  carrier  confinement  and 
the  lattice  potentials.  However,  they  are  relatively  inflexible  in  studying  non-linear  dynamical  properties 
in  the  presence  of  strong  dissipation,  such  as  is  present  in  the  electron-phonon  interaction  at  high  electric 
fields. 

Another  approach  to  quantum  transport  relies  on  the  "integral  over  paths"  method,  originally  pro¬ 
posed  by  Dirac  and  formulated  by  Feynman  [5].  Path  integral  methods  rely  on  an  influence  functional 
technique  in  which  the  source  of  the  dissipation  has  been  integrated  over  all  phonon  modes.  This  results 
in  a  model  influence  functional  where  the  phonon-scattering  dissipation  can  be  represented  as  an  interac¬ 
tion  with  a  collection  of  harmonic  oscillator  modes  in  which  the  translational  invariance  of  the  carriers  is 
preserved.  The  resulting  model  includes  constant  or  oscillatory  electric  and  magnetic  fields,  carrier 
screening,  scattering  and  dissipation,  carrier  confinement,  background  temperature  and  initial  conditions 
can  be  dealt  with  as  readily  as  for  a  free  particle  [6]. 

The  work  at  N.  C.  State  University  supported  by  the  ONR  has  progressed  over  the  past  fifteen  years 
from  the  realm  of  particle  models  to  quasi-particle  models  to  quantum  transport  models.  We  have  relied 
primarily  on  the  Monte  Carlo  method  to  study  and  solve  transport  problems  in  III-V  compound  semicon¬ 
ductors.  Our  techniques  for  modeling  materials  physics  and  device  phenomena  extend  to  device  dimen¬ 
sions  around  one  thousand  angstroms  [7-8],  and  we  are  confident  in  these  models  for  predicting  steady 
state  and  transient  device  effects  down  to  these  dimensions.  However,  in  order  to  remain  in  the  forefront 
of  transport  physics  and  device  research,  we  must  continue  and  increase  our  progress  into  the  realm  of 
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dimensions  where  physical  effects  can  be  studied  in  device  regions  with  dimensions  less  than  one 
thousand  angstroms. 

As  one  new  direction  for  our  research,  we  have  chosen  to  pursue  the  path  integral  method.  This  is 
in  contrast  to  other  methods  now  being  studied  and  supported  by  ONR,  such  as  the  application  of  the  den¬ 
sity  matrix  formalism  and  Wigner  distribution  function  approaches  including  the  use  of  moment  equa¬ 
tions  [9].  We  have  made  substantial  progress  in  the  path  integral  approach  during  the  past  four  years  and 
we  believe  that  path  integral  methods  will  play  an  increasingly  more  significant  role  in  developing  an 
understanding  of  quantum  transport  in  devices.  Currently,  we  are  aware  of  only  one  other  U.S.  univer¬ 
sity  program  directed  toward  applying  the  path  integral  method  to  semiconductor  devices.  This  is  the  pro¬ 
gram  at  the  University  of  Illinois  under  the  direction  of  Prof.  Karl  Hess.  Our  program  is  very  complimen¬ 
tary  to  the  one  at  Illinois,  and  we  are  in  fact  collaborating  with  those  personnel  through  the  use  of  the 
NSF  Supcrcompufing  Center  and  the  National  Center  for  Computational  Electronics  there.  One  of  ouf 
former  students  who  recently  finished  his  Ph.D.  degree  joined  Prof  Hess’  group  in  January,  1990  and  thus 
we  expect  this  interaction  to  increase.  We  propose  to  spend  a  major  effort  collaborating  with  Prof.  Hess’ 
group  in  applying  the  path  integral  to  specific  device  structures  such  as  mesoscopic  devices  [10]  during 
the  next  year,  thereby  further  establishing  its  utility  as  a  method  for  modeling  quantum  transport  in  dev¬ 
ices.  In  addition,  wc  will  continue  to  explore  the-Monte  Carlo  transport  method  for  study  of  novel  device 
structures  where  theoretical  underpinning  is  required.  One  particular  class  of  devices  which  we  will  con¬ 
tinue  to  study  is  based  on  hot  electron  injection  across  potential  barriers  (e.g.  heterojunctions  and  planar- 
doped  barriers)  where  the  transport  mechanisms  which  determine  device  behavior  arc  poorly  or  incom¬ 
pletely  understood.  Of  great  importance  here  is  the  fact  that  these  devices  offer  the  potential  to  compare 
experiment  and  theory.  Also,  these  devices  can  be  used  to  aid  in  our  understanding  of  quasi-ballistic  tran¬ 
sport  (3).  For  example,  we  have  explored  the  concept  that  ballistic-like  behavior  can  be  predicted  even 
when  .scattering  effects  arc  significant  [llj.  The  type  of  scattering  and  the  range  of  carrier  energies  are 
critical  in  this  regime.  Another  novel  class  of  structures  which  will  also  be  studied  in  detail  by  the  Monte 
Carlo  method  is  based  on  the  rapidly-emerging  pseudomorphic  or  strained-layer  devices.  These 
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structures  permit  extended  compositional  ranges  and,  thus,  have  a  number  of  potential  advantages  such  as 
higher  transconductance  and  channel  carrier  density.  Recent  studies  on  strain-induced  piezoelectric  fields 
have  opened  yet  another  possibility  in  realizing  ultra-fast  switching  devices,  resulting  in  an  increased 
importance  on  the  study  of  carrier  transport  in  these  devices.  Other  novel  devices  which  are  currently 
being  analyzed  include  delta-doped  field-effect  transistors  and  quantiun  well  structures,  based  on  studies 
which  have  begun  during  the  past  year.  Finally,  we  have  initiated  research  which  incorporates  Monte 
Carlo  methods  into  the  quasi-particle  approach  based  on  moments  of  the  Boltzmann  transport  equation 
with  quantiun  mechimical  corrections.  Preliminary  results  show  the  utility  of  this  approach  for  modeling 
microwave  and  millimeter  wave  devices  which  are  important  to  the  DoD’s  MMIC  program. 

We  want  to  re-iterate  the  important  interactions  which  have  developed  between  our  research  group 
and  other  researchers  during  the  last  year.  These  interactions  increase  the  impact  of  this  program  on  the 
field  of  semiconductor  device  physics  and  allow  our  research  efforts  to  be  far  more  productive  through 
increased  intellectual  efforts  and  enhanced  facilities  and  resources.  During  the  1990  funding  period,  we 
have  intensified  our  research  collaborations  with  Prof.  K.  Hess  at  the  University  of  Illinois,  Dr.  H.  Grubin 
of  Scientific  Research  Associates,  and  Drs.  M.  Stroscio  and  G,  lafrate  of  the  U.  S.  Army  Research  Office. 
These  research  collaborations  have  resulted  in  several  joint  publications  based  on  mutual  research 
interests,  expertise  and  capabilities.  There  have  been  several  visits  between  these  laboratories  and  the 
logistics  for  increased  collaborations  are  excellent. 

2.2  Summary  of  Research  Results 

2.2.1  Quantum  Transport 

This  section  will  be  divided  into  two  parts.  The  first  section  will  provide  a  summary  of  the  current 
status  of  our  work  on  the  development  of  Path-Integral  Monte  Carlo  (PIMC)  methods  for  the  study  of  the 
electronic  properties  of  ultrasmall  devices.  The  application  of  PIMC  is  based  on  the  Feynman  path- 
integral  (FPI)  formalism  of  quantum  mechanics.  The  second  section  will  provide  a  demonstration  of 
PIMC  methods  in  the  calculation  of  carrier-phonon  coupling  effects.  In  particular,  carrier  self-energies 
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are  calculated  by  PIMC  techniques  for  a  single-crystal  semiconductor  (GaAs)  and  for  GaAs  quantum 
wires  as  a  function  of  wire  size. 

A.  Summary  of  Path-Integral  Monte  Carlo  (PIMC)  Research  For  Ultra-Small  Device  Applications 

A.I.  Introduction 

Proper  treatment  of  scattering  processes  such  as  charge-carrier-phonon  coupling  is  crucial  for  accu¬ 
rate  modeling  of  the  electronic  properties  of  semiconductor  devices.  A  powerful  and  conceptually 
appealing  treatment  of  carrier-phonon  coupling  is  provided  by  the  Feynman  Path-Integral  (FPI)  formal¬ 
ism  of  quantum  mechanics  [5,12].  The  power  of  the  FPI  formalism  in  treating  coupling  processes  was 
demonstrated  tirst  in  quantum  electrodynamics  (QED)  [13],  but  soon  its  advantages  in  the  formally  simi¬ 
lar  treatment  of  carrier-phonon  coupling  in  bulk  semiconductors  also  was  realized  [14].  Further,  the 
reduction  in  device  sizes  to  dimensions  on  the  order  of  the  thermal  de  Broglie  wavelength  of  charge  car- 
tiers  offers  few  conceptual  difficulties  for  the  full  quantum  mechanical  FPI  formalism,  while  greatly 
reducing  the  value  of  the  seraiclassical  treatment  of  carrier-phonon  scattering  as  spatially  localized, 
instantaneous  events.  However,  the  reduction  of  device  dimensions  does  pose  great  practical  problems  to 
application  of  the  FPI  formalism.  For  analytical  evaluation,  both  the  device  potentials  and  the  net  effects 
of  coupling  potentials  must  be  modeled  as  either  linear  or  quadratic  in  the  carrier  coordinates  to  make  the 
calculations  tractable.  Thus,  while  a  variational  technique  introduced  by  Feynman  [14]  has  been  used 
with  great  success  to  model  carrier-phonon  coupling  in  bulk  semiconductors  [15-20],  extension  of  the  FPI 
formalism,  analytically,  to  the  study  of  ultrasmall  stmetures  has  been  slowed,  though  not  halted  [21-23], 
by  approximations  required  to  model  built-in  device  potentials. 

The  alternative  to  analytical  evaluation  of  FPIs  is,  of  course,  numerical  evaluation.  In  this  section, 
the  feasibility  of  numerical  evaluation  of  FPIs  in  ultra-small  device  applications  is  considered.  In  Section 
A.2  the  FPI  formalism  is  reviewed  briefly.  In  Section  A.3  a  practical  path-integral  Monte  Carlo  (PIMC) 
method  for  equilibrium  calculations  is  described.  In  Section  A.4,  various  basic  nonequilibrium  PIMC 
methods  and  their  limitations  are  presented.  Because,  to  the  best  of  our  knowledge,  no  nonequilibrium 
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PIMC  method  has  proven  widely  applicable  to  ultrasmall  device  modeling,  the  goal  in  this  latter  section 
is  to  exhibit  the  fundamental  problems  yet  to  be  overcome. 

A.2.  FPI  Formalism 

In  quantum  mechanics,  the  time  evolution  of  a  particle’s  wave  function  4'(x,t)  can  be  obtained  from 

+eo 

tb)  =  I  K(Xb,tb:x3,ta)'P(Xa,ta)dXa  where  K  (Xb,tb:Xa,ta)  is  the  quantum  mechanical  propagator  or 
Green’s  function.  In  the  FPI  formalism  K  (Xb,tb:Xa,ta)  is  given  by 


K(Xb,tb:Xa,la)  =  /expj  j  (-jmx^ft)  -  V[x(t)])  dtl  Dx(t) 


Xb 

where  — mx"(t)  is  the  kinetic  energy,  V[x(t)]  is  the  potential  energy  and  the  "  J  • '  •  Dx(t)"  is  a  notationai 

convenience  signifying  an  integral  over  all  conceivable  paths  x(t)  from  point  x,  to  point  Xb  [5,12],  To 
evaluate  Eq.  (1),  the  integral  over  paths,  the  "path-integral",  must  be  subdivided;  i.e.,  the  paths  must  be 
broken  down  into,  for  example,  straight  line  segments  or  Fourier  components  [5],  Also,  by  considering 
imaginary-time  arguments  in  the  propagator,  it  can  be  shown  that  the  equilibrium  density  matrix  subject 
to  a  Maxwell-Boltzmann  distribution  in  energy  at  temperature  T  is  given  by 


Xb 

p(Xb,Xa;T)  =  jexp- 

X. 


1  1 

-J  /  (ymx'-t-V[x(t)])dt[Dx(t) 


(2) 


where  kg  is  Boltzmann’s  constant  [5].  For  both  equilibrium  and  nonequilibrium  calculations,  the  net 
effects  of  coupling  charge  carriers  to  harmonic  systems,  e.g.  phonons,  can  be  incorporated  exactly  into 
path-dependent  influence  functionals  inside  the  path-integrals  [5],  Though  this  properly  is  largely  the 
motivating  force  for  developing  numerical  FPI  tecluiiques  for  ultrasmall  device  applications,  inclusion  of 
influence  functionals  is  not  a  primary  obstacle  to  the  development  of  these  techniques  and  their  considera¬ 
tion  here  would  complicate  the  discussion  needlessly. 


A.3.  Equilibrium  (Imaginary-Time)  PIMC 


For  numerical  evaluation  of  FPIs,  to  achieve  a  reasonable  degree  of  resolution  in  the  path 
configiuations,  the  number  of  path  segments  or  Fourier  components  used  to  describe  the  paths,  though 
finite,  must  be  large.  The  resulting  high  dimensional  (path)  integral  inherently  is  no  amenable  to  evalua¬ 
tion  by  uniform  sampling  techniques  and,  thus,  path-integral  Monte  Carlo  (PIMC)  methods  must  be  used. 
Though  there  are  a  variety  of  methods  for  numerically  evaluating  imaginary-time  FPIs  [24],  the  method 
most  amenable  to  ultrasmall  device  applications  appears  to  be  direct  sampling  [25,26].  Paths  Xn(t)  are 
selected  stochastically  with  a  probability 


p[Xn(t)]  exp^ 


1  1 

“¥  J  ymin(0dtk 


(3) 


ft 

The  equilibrium  density  matrbc  then  is  obtained  from 


1  N 

p(Xb.x„:T)  =  pf(Xb,Xa:T)—  Xexp^ 


K/IcbT 

I  V[x„(t)]dt^ 


(4) 


where  P(<Xb,x.,;T)  is  the  free-space  density  matrix  and  N  is  the  total  number  of  sample  paths.  Despite  -  or, 
more  accurately,  because  of  -  the  simplicity  of  this  method,  it  has  several  properties  that  allow  its  practi¬ 
cal  application  to  modeling  equilibrium  charge  carrier  behavior  in  ultrasmall  devices;  (1)  Either  in  terms 
of  Fourier  coefficients  of  using  a  sequential  bisection  technique  [27],  paths  can  be  selected  from  Gaussian 
probability  distributions  with  no  unwanted  correlation  between  successively-generated  paths.  (2)  Paths 
can  be  stored  and  reused  to  model  various  potential  structures  and  can  be  normalized  readily  with  respect 
to  temperature,  and.  thus,  the  generation  of  paths  is  in  effect  a  start-up  cost.  (3)  By  reusing  paths,  small 
changes  in  the  density  matrix  between  similar  potential  structures  can  be  evaluated  accurately,  including 
changes  much  smaller  than  the  statistical  fluctuations  obtained  in  evaluation  of  each  structure  indepen¬ 
dently.  (4)  Energy  values  can  be  obtained  directly  from  the  probability  density  normalization  constants. 
(5)  Because  it  is  a  Monte  Carlo  technique,  little  additional  cost  is  incurred  in  treating  potentials  that  vary 
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in  multiple  dimensions,  e.g.  the  charge-carrier-phonon  coupling  potential.  With  this  equilibrium  PIMC 
method,  for  example,  energy  level  splittings  in  coupled  quantum  well  systems  of  a  few  percentage  points 
of  the  nominal  value  can  be  accurately  resolved  in  about  an  hour  on  a  VAXstation  3200  [26].  Though 
such  calculations  can  be  perfoimed  more  readily  using  other,  non-PIMC  techniques,  these  latter  tech¬ 
niques  are  not  as  amenable  to  physically  accurate  treatment  of  charge-carrier-phonon  coupling  processes 
in  ultrasmall  devices.  With  PIMC,  however,  it  should  be  a  straightforward  procedure  to  introduce  these 
coupling  processes  into  future  calculations. 


A.4.  Nonequilibrium  (Real-Time)  PIMC 

Though  nonequilibrium  and  equilibrium  FPI  calculations  pose  similar  analytical  computational 
problems,  nonequilibrium  PIMC  calculations  pose  far  more  severe  numerical  problems  than  do  equili¬ 
brium  PIMC  calculations.  For  nonequilibrium  calculations,  because  the  argiunent  of  the  exponential  in 
Eq.  (1)  is  purely  imaginary,  these  is  no  inherent  normalizable  probability  distribution  from  which  paths 
can  be  chosen  stochastically.  A  variety  of  methods  have  been  tried  to  overcome  this  fundamental  prob¬ 
lem.  The  most  direct  is  to  impose  a  decaying  probability  distribution  for  selecting  paths  away  from  clas¬ 
sical  paths,  and  then  increase  the  spread  in  the  distribution  function  until  the  solution  converges  [28], 
Another  functionally  related  method  is  suggested  by  rewriting  Eq.  (1)  in  the  equivalent  form, 


Xl,  Xb 


K(Xh,tb:Xa.t.O  =  11  w[x'(t)  -  x(t)]  exp^ 


•  1 

^J(|mx'2(t)-V[x'(t)])dt 


D/(t)Dx(t), 


(5) 


where  w[x'(t)-x(t)]  is  a  normalized  window  in  x'(t)  about  x(t)  (29,30).  The  window  is  taken  small 
enough  that  the  potential  function  may  be  approximated  either  linearly  or  quadratically  within  the  win¬ 
dow.  Therefore,  the  path-integral  over  x'(t)  can  be  performed  analytically.  The  result  of  the  first  integral 
is  a  function  of  x(t)  that,  not  coincidentally,  decays  away  from  classical  paths  providing  a  normalizable 
probability  distribution  for  selecting  paths.  Therefore,  the  second  integral  can  be  performed  using  PIMC 
techniques.  With  both  of  these  methods  the  absolute  value  of  the  contribution  for  each  sample  path,  an 
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exponential  in  a  purely  imaginary  argument,  is  independent  of  the  resolution  in  the  path  configurations, 
i.e.  of  the  number  of  path  segments  or  Fourier  terms  per  path.  However,  also  with  both  methods,  the 
average  value  of  these  contributions  decreases  geometrically  with  linear  increases  in  the  resolution  in  the 
path  configurations.  Consequently,  the  required  number  of  Monte  Carlo  samples  increases  geometrically 
with  linear  increases  in  path  resolution,  limiting  these  methods  to  short  time  transient  applications.  A 
fundamentally  different  approach  is  provided  by  viewing  real-time  FPIs  as  the  analytic  continuation  of, 
for  example  but  not  restricted  to,  imaginary-time  FPIs.  For  one  such  method  [31,32],  first,  values  of  a 
FPI  along  the  negative  imaginaiy-time  axis,  where  the  FPI  takes  the  form  of  Eq.  (2)  but  with  B/hgT 
replaced  by  the  magnitude  of  the  imaginary-time  period,  are  found  using  equilibrium  PIMC  methods. 
Next,  some  type  of  curve  fit  to  the  imaginary-time  solutions  is  performed.  Last,  an  analytic  continuation 
back  to  the  real-time  axis  is  made  for  the  curve  fit.  However,  this  approach  has  been  limited  by  instabili¬ 
ties  that  arize  during  the  analytic  continuation  of  the  curve  fit,  though  not  necessarily  from  the  curve  fit  a# 
evidenced  by  the  next  method.  The  next  and  final  nonequilibrium  PIMC  method  presented,  another  ana¬ 
lytic  continuation  method  but  one  that  requires  no  curve  fit,  is  analogous  to  and  can  be  derived  from  the 
direct  sampling  equilibrium  PIMC  method  described  in  .Section  A.3.  Path  variations  x'„(t)  from  the  con¬ 
stant  velocity  path  x„(t)  from  Xj  to  xj,  are  chosen  with  a  probability 


p[x'(t)]  oc  exd 


_1_ 

K 


The  real-time  FPI  then  is  obtained  from 


1  N 

K(Xb.tb;Xa.ta)  =  KKXb.tb:X.vta)—  X 


-^Jv[Xo(t)-i-x;(t)e'''/2]dtl 


(6) 


(7) 


where  KjfXb.tblXa.tj)  is  the  free-space  propagator  [33,34].  Here,  analytic  continuations  are  performed  for 
each  sample  before  summation.  Beyond  returning  the  time  argument  to  the  real-time  axis,  the  effect  of 
the  analytic  continuations  is  to  displace  the  previously  selected  path  variations  into  complex  position 
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space  where,  presumably,  the  analytic  continuation  of  V(x)  is  known.  The  value  of  the  real-time  FPI,  in 
theory,  then  can  be  found  by  summing  the  results  of  the  individual  analytic  continuations.  Unfortunately, 
for  potential  functions  that  are  analytic  eveiy  where  in  the  complex  plane,  as  the  size  of  the  path  variations 


x'„(t)  increases,  the  magnitudes  of  e.xp| 


K/ksT 


I  V[Xo(t)  +  Xn(t)e'’^]dt|-  tend  to  increase  much  faster  than 


the  probabilities  for  selecting  the  path  variations  decrease.  Therefore,  excluding  linear  and  some  qua¬ 
dratic  potentials  -  for  which  analytic  solutions  exist  already  -  and  near  diagonal  elements  of  the  propaga¬ 
tor  for  certain  barrier  potentials  for  short  time  periods,  the  summation  of  Eq.  (7)  does  not  converge  abso¬ 
lutely  -  and  possibly  not  at  all.  Thus,  again  the  numerical  solutions  are  unstable.  An  interesting  alterna¬ 
tive  would  be  to  use  potentials  functions  that  are  analytic  and  bounded  everywhere  in  the  complex  plane 
except  at  isolated  poles;  however,  it  is  not  clear  what  theoretical  and  practical  effects  the  resulting  intro- 
duction  of  essential  singularities  into  the  calculations  would  have. 


A.5.  Conclusion 

With  existing  methods,  PIMC  offers  a  practical  approach  to  modeling  the  equilibrium  electronic 
properties  of  ultrasmall  devices,  an  approach  that  should  provide  physically  accurate  treatments  of 
carrier-phonon  coupling  processes  in  the  future.  In  contrast  fundamental  computational  problems  so  far 
limit  existing  nonequilibrium  PIMC  methods  to  short-period  transient-time  applications.  While  future 
increases  in  computer  power  should  allow  nansient-time  applications  of  somewhat  longer  time  periods 
with  existing  methods,  it  appears  that  new  PIMC  methods  must  be  found  to  fully  overcome  these  prob¬ 
lems  and  allow  modeling  of  long-period  transient-tune  and  steady-state  nonequilibrium  behavior.  How¬ 
ever,  the  great  advantages  that  would  be  offered  by  a  practical,  widely  applicable  nonequilibrium  PIMC 
method,  for  modeling  charge-canicr  transport  in  ultrasmall  devices  as  well  as  in  other  areas  of  research. 


motivate  a  continued  search. 
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B.  Self-Energy  Calculations:  An  Example  of  the  Application  of  PMC  Methods 

As  a  demonstration  of  the  use  of  the  PIMC  method  to  calculate  carrier-phonon  coupling  effects,  in 
this  section  carrier  self-energies  are  calculated  for  free-space  (single  crystal)  as  a  function  of  temperature, 
and  for  "quantum  wires"  as  a  function  of  wire  size  at  77.3K  and  300K.  After  Degani  and  Hipolito  [34], 
the  self  energy  S  at  linite  temperatures  is  defined  by 

Z  =  <E>-<E^>  (8) 


where  <E>  and  <E„>  are  the  polaron  and  uncoupled  carrier  expected  energies  as  defined  by 

“■"A'*' 

^IcbT^ 


(9) 


where  F  is  the  free  energy  of  the  system  [5],  kg  is  Boltzmann’s  constant,  and  T  is  the  absolute  tempera- 
ture.  The  values  of  <E>  and  <Eo>  are  calculated  from  Eq.  (9)  using  Metropolis  importance  sampling  and 
a  finite  difference  approximation  [26]. 

Figure  1  shows  the  PIMC  calculated  free-space  self-energies  of  electrons  in  GaAs  as  a  function  of 
temperature  from  lOK  to  lOOOK.  The  solid  curve  is  for  only  one  path  set,  but  the  error  bars  at  77.3  and 
300K  show  the  small  (0.12%  and  0.24%  respectively)  RMS  deviation  based  on  10  path  sets  at  these 
points.  Because  there  was  no  varying  potential  energy  function  to  consider,  the  epu  time  required  for 
each  such  calculation  beyond  that  required  to  calculate  the  influence  functionals  was  negligible.  The 
dashed  curve  at  low  temperatures  keT  <  hto,  shows  results  based  on  the  variational  method  calculations 
of  finite  temperature  free-cnergies  of  Saitoh  [36].  The  calculations  appear  to  converge  at  very  low  tem¬ 
perature.  but  the  PIMC  results  demonstrate  a  .somewhat  more  rapid  falloff  in  self-energy  with  increasing 
temperature  that  the  variational  results,  though  in  either  case  the  falloff  Is  extremely  slow.  The  qualita¬ 
tive  structure  of  this  self-energy  vs.  tetnperature  curve  has  been  exhibited  previously  by  Degani  and 
Hipolito  [35]  using  a  variational  FPI  method;  however,  quantitative  comparison  is  not  possible  because 
their  results  were  given  for  quantum  wire  structures,  to  be  considered  below. 


SELF-ENERGY  (/iw) 


TEMPERATURE  (K) 


Figure  1.  Electron  free-space  self  energies  as  a  function  of 
temperature.  The  solid  line  represents  PIMC  results  for  a  set  of 
1000  paths.  The  dashed  line  represents  results  based  on  the 
variational  calculations  of  low-temperature  free-energy  by 
Saitoh  [36]. 
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Figure  2a  and  2b  display  the  means  and  RMS  deviations  in  PIMC  calculated  self-energies  of  square 
quantum  wire  structiues  (carrier  confinement  in  two  dimensions)  of  widths  varying  from  OA  and  200A. 
The  quantum  well  depth  was  taken  to  be  0.2158eV.  Each  such  calculation,  a  combination  of  two  separate 
free-energy  calculations,  required  approximately  two  hours  of  cpu  time  on  the  VAXstation  3200  at  77.3K 
and  one-half  hour  at  300K.  For  reference,  the  isolated  closed  circles  in  Figs.  5.2a  and  5.2b  are  OK, 
infinite  well  depth  results  obtained  analytically  by  Degani  and  Hipolito  [37]  using  a  variational  non-FPI 
method.  The  large  RMS  deviations  in  the  PIMC  results  at  77.3  for  the  midrange  of  well  widths  is  the 
result  a  large  majority  contribution  to  the  probability  density  from  a  small  minority  of  the  stochastically 
sampled  paths.  However,  again,  the  RMS  deviations  shown  here  are  for  sets  of  only  1,000  paths  each  and 
can  be  reduced  readily  by  increasing  the  number  of  paths  used;  the  root  of  the  expected  mean-square 
deviation  of  the  man  found  here  for  all  ten  path  sets  from  the  true  mean  is  one-third  that  shown.  At  ftrat 
for  both  temperatures,  judged  by  the  means,  the  PIMC  obtained  self-energies  increase  as  the  well  width 
decreases  and  the  electrons  become  more  localized.  However,  as  the  well  narrows  further,  the  self- 
energies  begin  to  decrease  again,  the  result  of  reduced  electron  confinement  as  the  electron  energy  levels 
are  forced  up  and  out  of  the  well.  The  oA  limit  is  just  the  free' space  condition  again  as  the  quantum  well 
vanishes.  This  ultimate  decrease  in  self-energy  with  well  width,  of  course,  is  not  obtained  under  an 
infinite  well  depth  approximation  for  which  carrier  confinement  to  the  well  is  complete  regardless  of  the 
wire  dimensions  [35.371.  As  can  be  seen,  judged  by  the  means  and  allowing  for  the  0.0034  Ii(o  and' 
0.0075  K(o  free-space  energy  shifts  from  OK  to  77.3K  and  300K,  respectively,  exhibited  in  Fig.  1,  the  iso¬ 
lated  OK  non-FPI-method  results  of  Degani  and  Hipolito  agree  well  with  the  PIMC  results  obtained  here. 
However,  the  agreement  of  their  variational  FPI-method  results  mentioned  above  to  both  the  PIMC 
results  and  their  own  non-FPI  results  is  poor.  In  particular,  inexplicably,  at  a  width  of  lOOA  their  FPI- 
methbd  results  at  both  OK  and  77.3K  exhibited  no  significant  confinement  effects  on  the  electron  self¬ 
energies,  in  marked  contrast  to  both  the  PIMC  results  and  their  own  non-FPI-method  results.  This  lack  of 
agreement,  perhaps,  can  be  considered  a  testament  to  the  difficulties  involved  in  applying  the  FPI  formal¬ 
ism  analytically  in  ultrasmall  devices  structures. 


SELF-ENERGY  (Aw) 


WIRE  WIDTH  (A)  WIRE  WIDTH  (A) 

(a)  (b) 


Figure  2.  Mean  and  RMS  deviations  of  PIMC  calculated  self¬ 
energies  in  square  quantum  wires  as  a  function  of  width  at  (a) 
77K  and  (b)  300K,  with  isolated  analytically-obtained  results  at 
OK  (open  circles)  after  Degani  and  Hipolito  [35]. 


15 


2.2.2  Monte  Carlo  Study  of  Ballistic  and  Quasi-Ballistic  Transport  in  Semicngductors 

This  section  summarizes  our  research  which  deals  with  electron  transport  across  potential  barriers. 
We  have  studied  two  particular  types  of  barriers,  which  are  the  planar-doped  barrier  and  the  abrupt 
heterojunction  barrier.  While  much  of  the  following  discussion  focuses  on  the  planar-doped  barrier,  the 
results  are  more  generic  and  seem  to  apply  to  potential  barriers  in  general,  at  least  where  the  potential 
changes  (variations)  occur  over  distances  smaller  than  the  electron  mean  free  path  at  the  total  barrier 
energy.  The  general  nature  and  applicability  of  the  physical  effects  and  resultant  physical  models  merit 
further  investigation. 

The  distribution  function  experiences  fundamental  changes  when  an  ensemble  of  carriers  passes 
through  an  abrupt  potential  barrier.  By  imposing  continuity  of  total  particle  energy  and  the  parallel  com¬ 
ponent  of  wave  vector,  there  are  two  important  consequences  on  the  velocity  distribution  function.  First, 

ft 

the  normal  component  of  the  velocity  vector  is  shifted  to  high  positive  values  by  an  amount  equal  to  the 
ballistic  velocity;  that  is,  the  velocity  gained  by  a  particle  accelerated  without  collisions  through  a  poten¬ 
tial  energy  equal  to  the  barrier  height.  Secondly,  the  width  of  the  distribution  is  reduced  by  the  non- 
parabolicity  of  the  energy  bands  and  by  the  reduction  of  the  effective  mass  in  crossing  the  barrier  (pro¬ 
vided  this  is  the  case).  This  effect  is  substantial,  producing  an  extreme  non-equilibrium  distribution  func¬ 
tion  with  excess  energy  (for  barrier  heights  well  above  the  thermal  energy)  and  a  large  forward  com¬ 
ponent  of  momentum.  This  situation  is  much  the  same  for  realistic  barriers.  For  abrupt  heterojunctions 
with  interfacial  space  charge  layers,  the  injected  velocity  distribution  function  is  essentially  indistinguish¬ 
able  from  an  ideal  abrupt  barrier.  For  planar-doped  barriers,  where  carrier  injection  occurs  over  a  finite 
high  field  region  a  few  hundred  angstroms  thick,  the  velocity  distribution  function  differs  only  slightly 
from  an  ideal  abrupt  barrier.  In  this  case,  there  is  a  "tail"  in  the  velocity  distribution  which  decreases 
toward  zero  from  the  ballistic  velocity.  Thus,  the  planar-doped  (or  delta-doped)  barrier  is  not  quite  as 
effective  in  transferring  the  barrier  energy  to  the  ensemble  due  to  scattering  events  in  the  high  field 
region.  However,  the  basic  nature  of  the  velocity  distribution  function  remains  phenomenologically  simi¬ 
lar  to  that  for  an  ideal  abrupt  barrier. 


The  particle  motion  away  from  these  barriers  is  not  controlled  by  normal  diffusive  processes  with 
diffusion  velocities  of  around  5  x  10®  cm/sec.  The  extreme  non-equilibrium  transport  allows  for  more 
than  an  order  of  magnitude  increase  in  the  ensem.ble  average  velocity  (~10*  cm/sec)  over  short  distances 
while  the  non-equilibrium  populations  are  reduced  to  related  (or  thermal)  distributions  through  isotropic 
scattering  events  which  tend  to  randomize  the  distribution  of  particle  velocities.  The  transport  of  the 
ensemble  during  this  transition  period  is  very  important  to  high-speed  device  behavior.  However,  it  is 
very  difficult  phenomena  to  incorporate  in  device  simulation  and  models. 

In  1985,  two  groups  reported  experimental  evidence  for  hot  electron  transport  in  barrier  injection 
devices  [38,39].  These  and  later  results  established  the  existence  of  a  peak  in  the  electron  energy  distri¬ 
bution  function  at  the  bjuxier  injection  energy  [40,41].  This  peak,  observable  a  significant  distance  from 
the  barrier,  was  attributed  to  ballistic  transport  in  both  homojunction  planar-doped  barriers  and  hctero- 
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Junction  barriers. 

Recently.  Hess  and  lafrate  explored  these  results  in  an  effort  to  elucidate  further  the  processes  of 
ballistic  transport  and  velocity  overshoot  [42].  They  developed  definitions  for  these  terms,  classified  the 
basic  experimental  results,  tmd  concluded  that,  while  the  basic  physics  of  ballistic  transport  is  well  under¬ 
stood  qualitatively,  much  quantitative  work  remains  to  be  done,  especially  with  respect  to  device  applica¬ 
tions. 

The  purpose  of  this  study  is  to  investigate  transport  in  the  base  region  of  GaAs  planar-doped  bar¬ 
riers  used  as  ballistic  launchers.  Oiu  results  expand  the  definitions  of  Hess  and  lafrate  and  further  classify 
electronics  in  the  velocity  distribution  function  as  to  their  role  in  ballistic  and  quasi-ballistic  transport  in 
semiconductors. 

The  basic  structure  of  a  planar-doped  barrier  has  been  discussed  at  length  in  the  literature  [3843]. 
Ideally,  the  device  has  a  triangular  potential  barrier  represented  by  two  regions  of  constant  electric  field 
strength.  A  plan-doped  barrier  transistor  (PDBT)  is  formed  by  constructing  two  such  barriers  separated 
by  a  base  region  with  zero  electric  field  strength.  From  Monte  Carlo  simulations,  the  electron  velocity 
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distribution  function  at  the  top  of  the  barrier  has  been  foimd  to  be  a  Maxwellian  distribution  with  essen¬ 
tially  no  negative  velocity  electrons  [43],  This  distribution  is  also  that  found  in  the  abrupt  emitter  of 
heterojunction  bipolar  transistors  (HBTs)  [44].  this  is  an  ensemble  effect  in  both  devices  since,  once  an 
electron  is  transported  over  the  barrier,  the  probability  for  back-scattering  is  very  small  [42-44],  These 
conditions  correspond  to  those  originally  described  by  Bethe  for  thermionic  current  in  Schottky  diodes 
[8],  and  have  been  verified  carefully  for  HBTs  using  Monte  Carlo  simulations  [44].  Thus  we  concluded 
that  the  emitter  current  of  a  PDBT,  as  well  as  an  HBT,  is  primarily  a  thermionic  current. 

The  velocity  distribution  function  at  the  edge  of  the  base  region  consists  of  two  components.  First, 
there  is  a  very  narrow  peak  centered  at  a  high  velocity.  This  velocity  equals  the  theoretical  free-flight 
velocity  calculated  for  an  electron  accelerated  in  a  non-parabolic  band  gaining  kinetic  energy  equal  to  the 
barrier  height  (field  time  distance).  Also,  there  is  a  small  tail  in  the  distribution  function  extending  from 
zero  velocity,  increasing  with  velocity  and  merging  with  the  narrow  peak.  This  is  very  similar  to  the  casS 
for  an  HBT  [44],  although  the  presence  of  the  tail  shows  that  the  planar-doped  barrier  is  not  as  effective  s 
the  abrupt  emitter  barrier  of  an  HBT  in  orienting  the  electron  velocity  in  the  transport  direction. 

As  the  electrons  iue  transported  from  the  emitter  into  the  base,  there  is  a  change  in  the  velocity  dis¬ 
tribution  function.  This  effect  is  shown  in  Figure  3  where  two  dominant  features  can  be  seen.  The  first  is 
the  peak  discussed  above  centered  at  the  free-flight  velocity.  This  narrow  peak  remains  for  a  substantial 
distance  into  the  base  beyond  the  mean  free  path  length.  However,  the  emergence  of  a  very  broad  peak 
centered  at  a  low  velocity  can  be  seen  also.  Thus,  there  is  a  redistribution  of  the  electrons  from  the  nar¬ 
row  peak  to  the  broad  peak,  suggesting  that  some  electrons  have  been  transported  far  into  the  base  with 
either  no  interactions  or  else  interactions  with  small  deviation  angle  and  energy  loss.  Other  electrons 
have  experienced  isotropic,  velocity-randomizing  interactions  leading  to  a  velocity-relaxed  population 
that  contributes  to  a  higlUy  diffusive  current  (Maxwellian  distribution).  Again,  the  same  phenomena 
occurs  in  an  HBT  [44],  Thus,  it  appears  that  electrons  in  the  base  of  a  PDBT  (and  an  HBT)  can  be 
divided  i.:to  three  generic  populations,  ballistic,  quasi-ballistic  and  velocity-relaxed.  This  classification  is 
consistent  with  the  definition  for  nearly-ballistic  transport  given  by  Hess  and  lafrate  [42]. 
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Figure  3.  Velocity  distribution  function  in  the  h^se  of  a  PDBT  at  a  distance  of  50  nm 
from  a  0.225  eV  barrier.  The  narrow  peak  centered  at  the  free  flight  velocity  was  the 
singular  velocity  distribution  at  the  edge  of  the  barrier;  the  broad  peak  was  not  present. 
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In  order  to  investigate  the  proposition  that  the  total  electron  distribution  can  be  divided  into  three 
populations,  we  have  performed  Monte  Carlo  simulations  using  the  following  three  definitions. 

i)  An  electron  is  in  the  ballistic  population  until  it  experiences  one  interactions. 

ii)  An  electron  is  in  the  quasi-ballistic  population  if  it  experiences  very  few  (^)  anisotropic  interac¬ 
tions.  As  anisotropic  interaction  is  defined  as  an  interaction  in  which  the  angle  between  the  initial 
and  final  electron  velocity  following  an  interaction  is  less  than  a  threshold  angle,  0,.  Here,  0,  and 
N  will  be  determined  empirically  from  the  Monte  Carlo  simulations. 

iii)  An  electron  is  in  the  velocity-relaxed  population  if  it  experiences  one  or  more  isotropic  interactions 
0  >  0a,  or  more  than  N  anisotropic  interactions. 

Monte  Carlo  simulations  with  several  values  of  0a  and  N  were  performed  and  the  velocity  distribution 
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functions  of  the  ballistic  (fj,  quasi-ballistic  (fj  to  fN)  and  velocity-relaxed  populations  (fj)  were  com¬ 
puted  50nm  into  the  base  region  from  a  0.225eV  barrier.  The  purpose  of  the  simulations  was  to  identify 
0a  and  N  such  that  the  narrow  peak  in  the  total  velocity  distribution  mainly  is  due  to  the  ballistic  and 
quasi-ballistic  populations,  and  the  broad  peak  is  due  to  the  velocity-relaxed  population.  Figine  4  shows 
the  results  of  these  simulations.  In  this  figure,  we  plot  the  area  under  the  function  [ft  -  (fo  +  ...  +  fn)] 
versus  N  for  different  values  of  0a,  where  f,  is  the  total  distribution.  Also  shown  is  the  area  under  a  Gaus¬ 
sian  function  lilted  to  the  broad  peak  in  the  total  velocity  distribution.  From  these  curves,  and  from  many 
other  simulations,  we  conclude  that  reasonable  empirical  values  are  N  =  3  and  0, «  70°.  The  results  of 
these  simulations  support  the  concept  that  it  is  possible  to  deconvolve  the  two  peaks  of  the  velocity  distri¬ 
bution  function  using  the  simple  criteria  described  above,  and  that  the  narrow  peak  is  indeed  due  to  ballis¬ 
tic  and  quasi-ballistic  electrons  while  the  broad  Gaussian  peak  is  due  to  velocity-relaxed  electrons.  It  is 
interesting  to  note  that  similar  simulations  of  HBTs  have  shown  that  the  contributions  to  the  base  ciurent 
from  the  ballistic  and  quasi-ballistic  populations  are  negligible  for  i  >  3  [44],  This  lends  further  support 
to  the  generic  definitions  for  the  ballistic  and  quasi-ballistic  populations  in  ballistic  launcher  structures. 
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and  extends  the  definitions  of  Hess  and  lafrate  [42]. 

An  analytical  model,  based  on  the  definitions  for  ballistic,  quasi-ballistic  and  velocity-relaxed  popu¬ 
lations,  has  been  developed  in  order  to  obtain  quantitatively  the  spatial  variation  of  the  various  physical 
variables.  The  model  is  based  on  the  decay  of  the  electron  flux,  Oi(x),  of  the  ith  ballistic  or  quasi-ballistic 
population  due  to  interactions  in  the  base  region.  By  developing  a  difference  equation  of  Oj,  we  can 
describe  the  exchange  of  flux  within  the  ballistic  and  quasi-ballistic  populations  for  each  differential 
region  of  the  velocity  distribution.  The  resulting  recursive  differential  equation  is 


dOjCz.Vj) 

dz 


- Oi(z,Vi)  +  —  <Di_,(z,Vi_,) 

V|  Vi_, 


(10) 


where  X,  is  the  total  scattering  rate  is  the  anisotropic  scattering  rate  and  vj  is  the  electron  velocity  of 

the  ith  population.  The  equation  for  each  value  of  i  can  be  solved  analytically  if  and  X„  are  assumed  to 
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be  constant.  Tltis  assumption  is  valid  for  the  energy  range  considered  here.  Thus,  if  we  know  the  value 
fo  the  flux  of  ballistic  electrons  z  =  0,  <l>o  (0,Vo),  we  can  solve  for  4>j(z,Vi)  for  i  =  1  to  3.  We  obtain 
^o(O.Vo)  from  the  expression  for  the  velocity  distribution  function  at  the  edge  of  the  base  [9],  Thus,  we 
can  construct  the  velocity  distribution  functions,  fj(z),  of  the  i  =  0  to  3  populations  for  all  values  of  z  in 
the  base.  The  electron  density,  ni(z),  average  velocity,  Vi(z),  and  current  density,  Ji(z),  of  the  i  =  0  to  3 
populations  are  calculated  from  the  velocity  distribution  functions.  The  total  ballistic  and  quasi-ballistic 
current  density,  J,,b(z).  is  then  given  by 

oo  3 

i=0  i=0 

Since  the  velocity  distribution  function  of  the  velocity-relaxed  electrons  is  Gaussian,  the  drift-diffusion 
formalism  can  be  used  to  calculate  the  population  density,  n,(z),  of  these  electrons.  This  analysis  results 
in  the  following  equations: 


dn,(z) 

Jr(z)  =  qDn  +  qPnEn,(z). 


(12) 
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Monte  Carlo  simulations  were  carried  out  to  validate  this  model.  Electrons  were  injected  from  the 
top  of  the  emitter  barrier  with  initial  momenta  selected  from  a  semi-Gaussian  function  (only  positive  z- 
directed  momenta).  An  ideal  collecting  plane  was  assumed  at  the  end  of  the  base  region.  Figure  5  shows 
the  close  agreement  obtained  between  the  model  and  Monte  Carlo  simulations  for  ni(z)Ji(z)  and  n,(z). 
The  values  of  the  total  scattering  rate,  Xf  and  the  anisotropic  scattering  rate,  Xa„,  obtained  from  these 
curve  fits  were  .1.18  x  lO'-’/sec  and  1.2  x  lO'^/sec,  respectively.  A  comparison  of  these  values  with  the 
total  scattering  rate  1.2  x  lO'Vsec)  used  in  the  Monte  Carlo  simulations  show  that  polar  optical  scatter¬ 
ing  is  indeed  the  primary  anisotropic  .scattering  process  in  the  range  of  energies  considered  here  (less  than 
-  0.3  eV),  and  that  the  model  can  be  used  for  any  material  once  the  total  and  anisotropic  scattering  rates 
are  identified  for  the  particular  study.  Figure  6  demonstrates  the  results  of  applying  such  a  model  to  a 
GaAs  planar  doped  barrier  with  a  0.225  eV  barrier  height.  This  figure  shows  that  the  current  flow  is  dom¬ 
inated  by  the  quasi-ballislic  transport  for  distances  much  greater  than  the  electron  mean  free  path  befoit 
the  traditional  drift-diffusion  current  begins  to  dominant.  This  same  effect  is  observed  from  Monte  Carlo 
simulations  of  heterojunctions  bipolar  transistors  [44].  In  general,  there  is  excellent  agreement  between 
the  model  and  Monte  Carlo  simulations. 

In  summary,  we  have  classified  electron  transport  in  the  base  of  a  planar-doped  barrier  transistor  in 
terms  of  ballistic,  quasi-ballistic  and  velocity-relaxed  populations.  Simle  criteria  have  been  developed  to 
define  these  populations  in  terms  of  the  degree  of  anisotropy  of  the  scattering  and  the  number  of  anisotro¬ 
pic  interactions  that  an  electron  experiences  during  transport.  Empirically,  an  interaction  is  classified  as 
anisotropic  if  the  electron  velocity  deviates  less  than  70°  due  to  the  interaction.  Also,  the  maximum 
number  of  anisotropic  interactions  which  an  electron  can  experience  and  remain  in  the  quasi-ballistic 
population  is  empirically  determined  to  be  equal  to  three.  These  definitions  allow  the  construction  of  a 
simple  imalytical  model  which  describes  the  spatial  variation  of  the  ballistic  and  quasi-ballistic  current  in 
the  base  region.  This  model  agrees  well  with  Monte  Carlo  simulations.  Despite  its  simplicity,  the  suc¬ 
cess  of  this  model  in  dc.scribing  the  decay  of  ballistic  and  quasi-ballistic  populations  makes  it  useful  in 
device  modeling  that  docs  not  require  intensive  computer  resources.  At  this  stage,  the  model  does  not 
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Figure  5a.  Comparison  of  analytical  model  with  Monte  Carlo  simulations  of  ballistic  0*o) 
and  quasi-bailistic  current  ’  densities  versus  distance  into  the  base  of  a  PDBT. 
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Figure  5c.  Comparison  of  analytical  model  with  Monte  Carlo  simulations  of  velocity - 
relaxed  electron  density  versus  distance  into  the^base  of  a  PDBT.  The  fitted  value  of  D; 
results  in  an  electron  mobility  in  the  base  of  8600  cm^/vsec. 
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Figure  5.  Current  versus  distance  in  a  GaAs  planar-doped 
barrier  with  a  barrier  height  of  0.225  eV.  The  data  points  are 
from  Monte  Carlo  simulations.  The  solid  curves  are  calculated 
from  the  analytical  model  described  in  the  text. 
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take  into  account  the  transfer  of  electrons  to  the  lateral  or  satellite  valleys,  Thus,  it  is  valid  only  for  the 
range  of  energies  where  intervalley  transfer  is  negligible.  However,  the  ideas  presented  here  can  be 
extended  to  include  the  effects  of  intervalley  scattering. 

2.2.3  A  Hydrodynamic  Hot-Electron  Transport  Model 

We  have  explored  a  new  approach  for  modeling  hot  electron  transport  in  submicron  III-V  com¬ 
pound  semiconductor  devices  during  the  past  three  years  which  we  believe  can  have  a  significant  impact 
on  the  development  of  device  models  for  microwave  and  millimeter  wave  circuit  simulators..  Currently, 
this  area  is  critical  for  the  DoD’s  MMIC  program,  and  there  is  an  urgent  need  for  computer-aided  design 
approaches  in  the  development  of  high  frequency  electronic  systems.  Often  the  design  of  such  circuits  is 
done  empirically  and  by  trial  and  error.  We  would  like  to  contribute  to  the  development  of  an  approach 
to  millimeter  wave  circuit  design  which  is  analogous  to  approaches  now  used  in  the  design  of  silicog 
integrated  circuits.  There  are  at  least  two  needs  for  the  implementation  of  such  an  approach.  First, 
models  are  required  for  the  description  of  high  frequency  circuit  effects  brought  about  by  the  interactive 
nature  of  circuits  elements  which  are  normally  considered  discrete  and  isolated  at  low  frequencies.  Thus, 
at  very  high  frequencies,  these  elements  cannot  be  considered  as  independent  circuit  elements  intercon¬ 
nected  by  lossless  wires,  but  their  mutual  interactions  must  be  considered  as  well  as  the  nature  and  physi¬ 
cal  layout  of  their  intercomiections.  Secondly,  specific  device  models  must  be  developed  which  take  into 
consideration  the  nonequilibrium  effects  which  occur  and  can  dominate  the  performance  of  submicron 
compound  semiconductor  devices.  Examples  are  quasi-ballistic  transport  and  velocity  overshoot,  which 
are  becoming  practical  limitations  of  current  device  structures  rather  than  the  esoteric  topics  they  were 
once  considered. 

We  arc  addressing  the  need  for  increased  physical  understanding  in  device  models  through 
advanced  modeling  concepts,  such  as  the  Feynman  path  integral  and,  on  a  more  practical  level,  through 
the  use  of  hydrodynamic  approaches  to  device  modeling.  This  (hydrodynamic)  approach  allows  none¬ 
quilibrium  effects  to  be  included  in  a  rigorous  way  within  device  models  which  are  more  readily  accepted 


by  circuit  and  device  engineers.  These  device  models  can  be  applied  for  dimensions  which  well  exceed 
the  current  technological  limits  of  device  fabrication. 

The  basic  transport  model  used  in  this  research  is  based  upon  the  seraiclassical  hydrodynamic  con¬ 
servation  equations  for  the  average  electron  density,  average  electron  momentum  and  average  electron 
energy.  The  general  model  includes  particle  relaxation  times,  momentum  relaxation  times  and  energy 
relaxation  times  as  a  function  of  average  carrier  energy  for  the  F,  L  and  X  conduction  band  valleys  of 
GaAs.  The  model  also  contains  electron  temperature  tensors  and  heat  flow  vectors  which  also  have  aver¬ 
age  energy  dependence.  The  relaxation  times  in  the  model  result  from  intravalley  and  interValley  particle 
.scattering  and  represent  rates  of  exchange  of  particles  between  valleys  and  rates  of  loss  of  average 
momentum  and  average  energy  between  and  within  the  individual  valleys.  The  electron  temperature  ten¬ 
sor  and  heat  flow  vector  depend  upon  the  electron  velocity  distribution  about  the  average  electron  velo- 

• 

city  and  ultimately  affect  transport  when  spatial  variations  in  average  velocity  and  average  energy  exist. 

In  the  initial  stages  of  the  project  [47],  the  extreme  heuristic  assumption  was  made  that  the  tempera¬ 
ture  tensor,  heat  flow  vector  imd  relaxation  times  were  all  simple  functions  of  average  energy.  These  gen¬ 
eral  transport  parameters  were  then  calculated  using  the  Monte  Carlo  method  in  conjunction  with  the 
ergodic  principle  applied  directly  to  the  integral  definition  for  the  parameters.  This  approach  was  taken  to 
approximately  include  the  effects  of  nonparabolicity  in  the  conduction  bands  and  to  make  as  few  assump¬ 
tions  as  possible  about  the  form  of  the  electron  distribution  in  momentum  space.  A  direct  or  indirect 
assumption  traditionally  made  is  to  assume  that  the  distribution  has  a  displaced  Maxwellian  form  which 
leads  to  incorrect  results  even  for  the  stationary  transport  case.  The  above  heuristic  approach  offered  a 
slight  improvement  in  accuracy  and,  if  developed,  offers  a  future  hydrodynamic/Monte  Carlo  model 
which  could  quickly  and  accurately  describe  many  electronic  structures. 

The  initial  investigations  of  the  hydrodynamic  model  have  been  based  on  single  conduction  band 
valley  cases  to  simplify  the  difficulties  associated  with  determining  solutions.  Solutions  were  obtained  for 
the  model  by  combining  a  very  efficient  local  Newton-based  solver  with  a  perturbation-in-doping  based 


23 


continuation  method.  Tfie  first  simulation  results  [48]  of  the  hydrodynamic  model  indicated  a  general 
improvement  over  the  more  conventional  drift-diffusion  model  with  results  closer  to  those  of  a  Monte 
Carlo  transport  model.  The  results  also  shewed  an  extreme  dependency  on  the  functional  form  of  the  heat 
flow  vector.  The  heat  flow  vector  was  shown  to  have  a  large  impact  on  the  resultant  cooling  in  average 
energy  observed  at  the  source  end  of  an  structure.  The  complete  results  from  this  study  are 

shown  in  Figs.  7-12  with  Fig.  12  indicating  an  excessive  cooling  at  the  source  which  did  not  agree  with 
the  Monte  Carlo  results.  Another  important  result  discovered  in  these  initial  investigations  was  an  inabil¬ 
ity,  due  to  neglecting  any  form  of  velocity  or  kinetic  energy  dependency  in  the  temperature  tensors,  to 
find  high  bias  solutions  for  the  model. 

Since  simulation  results  indicated  it  was  important  to  include  additional  accuracy  into  the  model 
through  more  physically  correct  effective  temperature  dependencies,  a  new  derivation  was  performed.  As 
in  the  previous  work  the  new  hydrodynamic  model  was  based  on  a  unique  set  of  moments.  The  usu^ 
approach  is  to  take  moments  of  the  electron  momentum  which  yield  a  set  of  conservation  equations  for 
density,  momentum  and  energy.  Then  addition  assumptions  are  used  to  mathematically  close  the  moment 
process  at  the  conservation  of  energy  equation.  The  traditional  assumption  is  to  either  directly  or 
indirectly  assume  a  displaced  Maxwellian  distribution.  Under  this  assumption  one  can  express  the 
moment  equations  in  terms  of  the  variables  density  n,  average  velocity  v,  average  energy  w,  and  thermal 
temperature  T.  The  displaced  Maxwellian  approximation  also  yields  the  constitutive  relation 

3  I  1 

w  =  resulting  in  a  value  of  the  heat  flow  vector  of  q  =  0  for  parabolic  bands.  Also,  a  heat 

flow  vector,  with  higher  order  differential  dependence  on  energy,  is  usually  included  using  thermodynam¬ 
ics  arguments  (phenomenological  Wiedemann-Franz  law). 

The  moments  used  in  the  new  model  can  be  described  as  moments  of  the  nonparabolic  electron 
velocity.  Specifically,  the  current  hydrodynamic  model  was  developed  using  the  moment  operators 

'Fo  =  l,  'F|  =  u(k)  and  'Ft  =  y  m(k)u(k)u(k)  =  Ec(k).  Here,  the  Kane  dispersion  relation, 

K“k“/2m*  =  Ec(l +aE<.)  is  used  to  define  the  single  electron  velocity  u(k)  and  the  effective  mass 
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Figure  11.  Drift  velocity 
profile  for  an 
applied  potential 
of  0.4  volts. 
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Figure  12.  Average  energy 
profile  for  an 
applied  potential 
of  0.4  volts. 
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m(k)  =  m*(l  -fc2aR‘kVin*)'^.  These  moment  operators  were  chosen  because  they  lead  to  a  form  which 
can  be  manipulated  more  easily,  one  in  which  simplifying  assumptions  can  be  seen  more  clearly.  Per¬ 
forming  the  moment  process  and  simplying  to  first  order  yields  immediately  the  results: 


=  -Vr-(nv) 


„  F  1  r7  r«  , 

—  =  -vV,V  +  — - rVr-[Pv] 

ot  m  nm 


=  -vV,w(PJ+q) 


where  m**,  [?v].  [Pwl  ‘>nd  q  have  integral  definitions  which  depend  on  f.  Here,  unlike  the  parabolic  band 
case,  there  is  a  dilemma  of  two  distinctly  different  effective  pressure  tensors  Py  and  P,,,  (for  parabolic 
bands  and  displaced  Maxwellian,  (PyJ  =  To  resolve  this,  stationary  Monte  Carlo  calculations  and  phys^ 
cal  intuition  arc  used  to  suggest  hie  following  ansatz.  distribution  function 


f  mw(T^v)  1  ■  ■  r  -mw(Tw)  I  u  -  v  I  ^ 

""  2nkBTw  2kBTvy 


as  a  constitutive  relation  to  close  the  moment  equations.  Here  T^  has  been  chosen  to  replace  T  because 

3 

— kfiTw  accurately  approximates  the  effective  thermal  energy  for  the  stationary  transport  case,  and 

hiw(Tw)  has  been  introduced  because  we  expect  a  non-constant  effective  mass  strongly  dependent  on  T^,. 
Equation  (16)  is  now  used  in  the  definitions  for  the  transport  parameters,  and  the  analysis  is  limited  to 
first  and  second  orders  in  T^,  to  arrive  at  the  the  supplemental  relations; 


w  =  ^kn'rw  + 


in"(TJ 


m“(TJ  =  nr(l+3akBTJ 


[Py)=~;;-nkBTJI] 

m 


(19) 


25 


[P^]  =  nkBTJI] 


qi  =  5a  —  (keTJ^nui 
m 


These  relations  supply  all  the  necessaiy  equations  to  express  the  hydrodynamic  transport  model  for  non¬ 
parabolic  conduction  bands  in  terms  of  the  quantities;  electron  concentration  n,  average  electron  velocity 
V,  average  electron  energy  w  and  the  space-times  variable  (r,t).  Applying  these  relations  reveal  the  new 
collisionless  nonparabolic  hydrodynamic  transport  model  equations  to  be : 


=  -Vr-(nv) 


d\  „  F  2  _ .  ,  m  . - 
—  =  -vV^v  +  ■—  -  - — rVr[nv(w  -  — vv)] 
m  |Li  3nni  ^ 


ft"  “  “  ■^^r'[nv(©  +  4>)(w  -  -^vv)] 


with  nonparabolic  correction  tenns 


li(v,w)  =  1  +  2a(0(v,w)(w  — 


v(v,w)  = 


0)(v.w) 

p(v.w) 


0)(v,w)  = 


1  +  otm*vv 


^  10a  oKv.w)  ^  m 

<^{v,w)  =  — - - — -  (w  -  vv) 

3  p(v.w)  2 


The  nonparaboiic  model  equations  above  have  been  written  in  a  form  so  that  for  the  parabolic  case,  a=0, 
it  is  easy  to  observe  that  p=v=a>=l  and  4)=0.  Then,  the  nonparabolic  model  conveniently  reduces  to  the 
more  familiar  classical  hydrodynamic  equations  in  Lagrangian  form. 
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Since  we  are  also  interested  in  models  which  can  be  used  to  study  ultrasmall  electronic  devices 
investigations  into  adding  quantum  corrections  to  the  hydrodynamic  equations  were  performed.  To 
achieve  this  goal,  the  previous  work  of  Grubin  and  Kveskovsky  [49]  was  followed  to  develop  a  set  of 
quantum  balance  equations  consistent  with  the  nonparabolic  model  above.  For  the  initial  investigation  a 
slightly  restricted  version  of  their  quantum  hydrodynamic  equations  were  used.  Specifically,  spatial  and 
nonparabolic  effects  on  the  effective  mass  were  neglected  and  any  deficiencies  of  the  general  quantum 
distribution  function  to  agree  with  Fermi  statistics  (low  temperature  effects)  were  ignored.  Development 
of  the  quantum  correction  terms  began  with  the  general  moment  equations  of  Stroscio  [50].  The  colli¬ 
sionless  one-dimensional  forms  of  the  first  three  moments  are: 


3n 

at 


JL± 
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(npd) 


(29) 


3pd  a  .  Pd  \ 


(30) 
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2Vx  P<i((P-Pd)^  an 


nm 
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(31) 


where  n  is  the  electron  density,  p  is  the  single  electron  momentum,  pj  is  the  classical  average  momentum  • 
and  Ueff  is  the  total  effective  electric  potential.  Here  Ref.  49  was  followed  and 


fw  =  — exp 


-P(P  -  Pd)’ 


2m 


I  Y  /I  P  /  x2  \  ^  1  \ 


(32) 


which  is  the  momentum  displaced  nonequilibrium  Wigner  distribution  function  of  Ancona  and  lafrate 

[51],  was  u.sed  to  obtain  n((p  -  pj)^  =  -^(1  -  "(^P  “  Pd)^= 0  where  a= li^P^/Sm*, 

P=  l/kgT  and  Y=cc/p.  Tliese  results  can  be  used  with  Eqs.  (29)-(31)  to  yield  the  collisionless  quantum 
hydrodynamic  equations: 
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where  w=  — IcbT  +  —  m  v- - r^(~-::“)  and  Q  = - -  r-  "•  Comparmg  the  above  quan- 

2  2  24m  9x '  n  9x '  ^  2m  ^  d\ 

turn  equations  to  the  semiclassical  nonparabolic  hydrodynamic  Eqs.  (22)-(24),  approximating  by  T 
and  using  the  supplemental  relations  (25)-(28)  one  can  arrive  at  the  quantum  corrected  collisionless  non¬ 
parabolic  hydrodynamic  transport  model  in  one  dimension, 
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with  nonparabolic  coixeclion  terms  p(n,v,w)  =  1  +  2a  co(v,w)  (w  ~  Wq - - — X  v(n,v,w)  =  < 
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]“(w-Wq - - — )  and  quantum  corrections  Wq 
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An  initial  study  has  been  performed  with  this  quantum  hydrodynamic  model.  The  quantxun  hydro- 
dynamic  model  was  used  to  analyze  a  double-barrier  heterostructure,  Al^pai.^^As  barriers  with  Xc=0.3, 
doped  to  10'‘*cm“^  with  source  and  drain  regions  of  length  200A  and  doped  to  10‘®cm“^.  This  structure, 
which  has  barrier  widths  of  sA  and  a  spacing  of  50A,  is  .shown  in  Fig.  14  with  the  source  and  drain 
regions  excluded.  In  the  self-consistent  simulations,  the  hetefostructure  was  treated  using  the  Anderson 
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treated  in  a  classical  manner  and  we  assumed  T|a„ice  =  300  K.  Simulation  results  were  generated  for  a 
number  of  biases  in  the  range  5  mV  to  50  mV  which  yielded  an  approximately  linear  J-'F^ppUed  charac¬ 
teristic  of  slope  0.36  kA/cm"-mV.  The  profile  result  of  w,  and  Wq,  parameters  intrinsic  to  the  energy  bal¬ 
ance  equation,  are  displayed  in  Fig.  13  for  'Fappiied=8mV.  Figure  14  shows  the  corresponding  profile  for 
Q  along  with  the  conduction  band  Ec  (Ec=Ueff). 

This  study  was  the  first  to  use  the  three  quantum  hydrodynamic  (QHD)  equations  to  study  transport 
in  heterostructurc  devices.  The  key  additional  feature  of  the  QHD  equations  is  the  incorporation  of  the 
effects  of  density  gradients  through  the  quantum  potential.  Insofar  as  a  classical  solution  does  not  exist  in 
the  presence  of  barriers,  Q  must  at  least  cancel  the  effects  of  the  barriers  and  permit  carrier  transport. 
Indeed  Q/3  does  approximately  balance  the  barrier  potential  a  result  consistent  with  Ref.  49.  Also,  since 
this  result  is  near  equilibrium,  the  density  should  vary  everywhere  as  r4*’exp[-(Ueff  +  Q/3)/kBT]  and  this  is 
found  to  be  the  case. 

2.2.4  Transport  in  Photoexcited  Semiconductors 

During  the  past  tluec  years,  we  have  begun  an  effort  to  study  the  transport  of  photoexcited  carriers 
in  .semiconductors.  We  have  decided  to  explore  the  use  of  Monte  Carlo  simulation  as  well  as  quantum 
mechanical  wave  approaches.  The  specific  device  which  we  have  chosen  to  first  explore  is  the  GaAs 
metal-.semiconducior-metal  (MSM)  photodetector  with  planar  (low  interelectrode  capacitance)  structure 
which  is  of  interest  for  discrete  applications  and  optoelectronic  integrated  circuits  (OEIC’s).  The  fastest 
reported  MSM  detector  has  an  intrinsic  bandwidth  of  105  GHz,  which  compares  well  with  btJidwidths  of 
other  Schottky  photodiodes.  These  low  noise,  high  quantum  efficiency  detectors  are  fully  MESFET  com¬ 
patible.  The  most  complex  high  performance  OEIC  uses  a  refractory-gate,  ion-implanted  process  with 
over  2000  devices.  Tlicre  is,  however,  little  simulation  work  done  for  devices  with  submicron  distance 
between  fingers.  Recently,  we  have  analyzed  detectors  with  dimensions  possible  to  obtain  with  well 
refined  lithography  techniques  [52-54].  Some  of  the  preliminaiy  results  on  temporal  response  and  dark 
current  were  reported  last  year.  In  this  section,  we  present  a  more  detailed  time  analysis  of  photodetec- 
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rule;  Ueff  =  q'Feff-x(x),  where  x(x)  (AE<.  =0.697Xc)  is  a  position  dependent  electron  affinity  and  the 
applied  potential  'F^ff  was  determined  from  solving  Poisson’s  equation.  The  dissipation  mechanisms  were 
treated  in  a  classical  manner  and  we  assumed  T|a„icc  =  300  K.  Simulation  results  were  generated  for  a 
number  of  biases  in  the  range  5  mV  to  50  mV  which  yielded  an  approximately  linear  J-SKappijed  charac¬ 
teristic  of  slope  0.36  kA/cm--mV.  The  profile  result  of  w,  and  Wq,  parameters  intrinsic  to  the  energy  bal¬ 
ance  equation,  are  displayed  in  Fig.  13  for  'Pappiied=8mV.  Figure  14  shows  the  corresponding  profile  for 
Q  along  with  the  conduction  band  E<.  (E<.=Ueff). 

This  study  was  the  first  to  use  the  three  quantum  hydrodynamic  (QHD)  equations  to  study  transport 
in  heterostructure  devices.  The  key  additional  feature  of  the  QHD  equations  is  the  incorporation  of  the 
effects  of  density  gradients  through  the  quantum  potential.  Insofar  as  a  classical  solution  does  not  exist  in 
the  presence  of  barriers,  Q  must  at  least  cancel  the  effects  of  the  barriers  and  permit  carrier  transport. 
Indeed  Q/3  does  approximately  balance  the  barrier  potential  a  result  consistent  with  Ref.  49.  Also,  since 
this  result  is  near  equilibrium,  the  density  should  vary  everywhere  as  I4’‘exp[-(Ueff + Q/3)/kBT]  and  this  is 
found  to  be  the  case. 

2.2.4  Transport  in  Photoe.xcited  Semiconductors 

During  the  past  tluee  years,  we  have  begun  an  effort  to  study  the  transport  of  photoexcited  carriers 
in  semiconductors.  We  have  decided  to  explore  (he  use  of  Monte  Carlo  simulation  as  well  as  quantum 
mechanical  wave  approaches.  The  specific  device  which  we  have  chosen  to  first  explore  is  the  GaAs 
metal-semiconductor-metal  (MSM)  photodetector  with  planar  (low  interelectrode  capacitance)  structure 
which  is  of  interest  for  discrete  applications  and  optoelectronic  integrated  circuits  (OEIC’s).  The  fastest 
reported  MSM  detector  has  an  intrinsic  bandwidth  of  105  GHz.  which  compares  well  with  bandwidths  of 
other  Schottky  photodiodes.  These  low  noise,  high  quantum  efficiency  detectors  are  fully  MESFET  com¬ 
patible.  The  most  complex  high  performance  OEIC  uses  a  refractory-gate,  ion-implanted  process  with 
over  2000  devices.  There  is,  however,  little  simulation  work  done  for  devices  with  submicron  distance 
between  fingers.  Recently,  we  have  analyzed  detectors  with  dimensions  possible  to  obtain  with  well 


Figure  13.  Average  energy  (w)  and  quantum  correction  (wq) 
for  an  applied  potential  of  8  mV. 
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Figure  14.  Conduction  band  edge  (Ec)  and  quantum  potential 
(Q)  for  an  applied  potential  of  8  mV. 
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refined  lithography  techniques  [52-541.  Some  of  the  preliminary  results  on  temporal  response  and  dark 
current  were  reported  last  year.  In  this  section,  we  present  a  more  detailed  time  analysis  of  photodetec¬ 
tors  with  0.1-1  pm  distance  between  fingers.  This  range  has  been  considered  because  it  spans  devices  of 
current  importance  as  well  as  those  of  possible  future  interest.  We  have  used  a  self-consistent  Monte 
Carlo  (MC)  model  to  simulate  intrinsic  transport  and  a  lumped-constant  circuit  model  for  parasitics.  The 
MC  results  are  used  as  an  input  for  simulation  of  the  circuit  model.  This  approach  is  relatively  easy  to 
implement.  However,  the  interpretation  of  the  results  depends  critically  on  the  parasitic  lumped-constant 
circuit  parameters. 

The  GaAs  MSM  photodetector  structure  has  been  illustrated  in  Ref.  53.  We  consider  two  different 
versions  of  this  detector:  di.scrcte  and  monolithically-integrated.  In  the  discrete  version,  the  bond  pads 
are  50  x  50  pm’  and  the  device  must  be  wire  bonded  to  the  load.  The  inclusion  of  the  series  wiring 
inductance  affects  detector  pulse  responses.  In  the  monolithic  version,  one  bond  pad  is  replaced  with  a 
microstrip  transmission  line  (10  x  50  pm’),  thus  reducing  the  series  inductance.  All  dimensions  are 
chosen  to  closely  model  realistic  devices.  A  low  bias  is  essential  if  the  photodetector  is  to  operate  in  a 
highly-integrated  optoelectronic  circuit.  It  is  fortunate  that  the  MSM  photodetector  offers  high  perfor¬ 
mance  at  relatively  low  bias  voltages.  Due  to  the  very  low  doping  concentration  of  the  epitaxial  layer,  the 
regions  between  and  below  fingers  are  completely  depleted  of  free  carriers  [53]. 

Photodetectors  with  separations  between  fingers  of  0.1-1  pm  were  simulated  and  results  are  shown 
in  Fig.  15.  A  self-consistent  ensemble  Monte  Carlo  simulation  was  used  to  separately  determine  electron 
and  hole  currents  for  a  large  number  (15000-25000)  of  photoexcited  carriers.  The  very  sharp  electron 
peak  for  the  device  with  a  0.5  pm  distance  between  lingers  is  due  to  a  higher  electric  field  than  in  other 
devices.  Although  the  intrinsic  re.spon.ses  were  calculated  for  the  detectors  with  a  lightly  doped  epitaxial 
layer,  one  can  expect  that  very  similar  results  would  be  obtained  if  the  devices  were  processed  directly  on 
a  semi-insulating  substrate.  We  observe  that  extremely  short  current  pulses  are  possible.  However,  other 
factors  have  a  strong  influence  on  output  signal  of  the  device. 
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Figure  15.  Comparison  of  intrinsic  responses  (no  parasitics) 
for  photodetectors  with  submicron  electrode  pattern  (<1  pm). 
As  interfinger  spacing  increases  from  0.5  pm  to  1  pm,  the 
longer  hole  tail  decays  become  more  evident. 
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The  parasitic  circuit  of  a  photodetector  (Fig.  16)  can  have  a  very  detrimental  effect  on  time 
response.  Briefly,  parasitic  elements  are  estimated  in  the  following  way.  Here,  CPD  is  modeled  as  a 
capacitance  between  parallel  microstriplines.  Also,  RF  and  LF  are  the  resistive  and  inductive  parts  of  the 
finger  impedance  (skin  depth  at  100  GHz  has  been  included).  In  the  discrete  version,  CPP  is  the  capaci¬ 
tance  between  bond  pads,  CPG  is  the  capacitance  between  the  bond  pad  and  ground,  and  LL  is  the  wire 
bond  inductance  (0.06  nH).  In  the  monolithic  version  LLl,  LL2,  CLG,  CLL  represent  the  microstripline. 
Other  parasitics  (RD.  CFG)  are  less  important,  but  were  included.  The  calculated  intrinsic  responses  can 
now  be  used  to  determine  the  responses  of  photodetectors  with  parasitics.  In  Fig.  17  we  present  the  nor¬ 
malized  output  voltages  vs  time  for  the  discrete  detectors.  The  FWHM  is  3.77  ps,  3.5  ps,  3.28  ps,  3.08  ps, 
2.86  ps  for  the  discrete  devices  with  spacings  between  fingers  of  0.1  pm,  0.25  pm,  0.5  pm,  0.75  pm,  and 
1  pm,  respectively.  Due  to  the  long  hole  tails  and  oscillations  present  in  some  devices  (monolithic),  the 
FWHM  can  not  be  used  unambiguously  to  estimate  the  bandwidth.  Therefore,  we  define  the  maximum 
repetition  frequency  as  follows.  If  t,ni„  is  the  time  for  the  pulse  to  reach  3%  of  maximum  (or  -3%  if  pulse 
becomes  negative),  and  t,„a^  is  the  time  to  reach  the  same  value  before  the  pulse  settles  to  zero,  then  the 
maximum  repetition  frequency  is  simply  l/ft^ax  “  imin)-  Based  on  this  definition,  the  maximmn  repetition 
frequency  is  calculated  and  results  are  shown  in  Fig.  18  for  the  intrinsic,  monolithic  and  discrete  detec¬ 
tors. 

The  experimental  results  related  to  the  MSM  photodetectors  are  somewhat  limited.  Although  there 
are  detectors  fabricated  with  a  submicron  electrode  pattern,  there  is  little  experimental  data  which  could 
be  used  for  the  purpose  of  comparison  'vith  simulated  results.  An  exception  to  this  situation  is  the  detec¬ 
tor  fabricated  and  tested  at  IBM  [551.  In  Fig.  19  we  present  the  experimental  data  and  the  simulated 
results  using  the  approach  described  above.  In  general,  there  is  good  agreement  between  these  two 
results.  However,  the  relatively  slow  leading  edge  of  the  experimental  pulse  is  difficult  to  explain.  This 
could  be  a  circuit  effect. 

Thus,  in  this  phase,  we  have  investigated  the  transit  behavior  of  the  GaAs  MSM  photodetectors. 
The  intrinsic  responses  have  been  simulated  using  Monte  Carlo  method  and  used  in  a  detailed  analysis  of 
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Figure  16.  Equivalent  circuit  of  a  MSM  photodetector.  The 
intrinsic  current,  IpD>  obtained  from  Monte  Carlo  simulation  of 
the  intrinsic  device,  is  used  to  excite  the  parasitic  circuit.  The 
parasitic  elements  depend  on  whether  the  discrete  or  monolithic 
version  of  the  photodetector  is  considered. 


Figure  17.  Time  responses  of  discrete  photodetectors  with 
parasitics.  The  full-width-half-maximum  (FWHM)  decreases  as 
interfinger  spacing  increases.  The  maximum  repetition 
frequency,  however,  is  determined  by  the  slow  hole  decay  tail 
rather  than  by  the  FWHM. 


max 


Distance  ({xm) 


Figure  18.  Maximum  repetition  frequency  versus  interfinger 
spacing  for  all  devices  simulated.  The  intrinsic  detector  is 
especially  fast  for  small  interfinger  spacing.  The  parasitic 
circuit)  however)  strongly  affects  the  maximum  repetition 
frequency. 
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the  discrete  and  monolithic  versions  of  the  detectors.  Instead  of  relying  on  the  FWHM  to  determine  the 
bandwidth,  we  defined  the  maximum  repetition  frequency  since  this  is  a  better  estimate  of  the  device 
speed.  Conducted  studies  reveal  that  the  GaAs  MSM  photodetectors  are  indeed  suitable  for  very  fast 
detection  of  optical  signals  in  both  di.screte  and  monolithic  versions.  These  studies  also  reveal  the  utility 
of  the  Monte  Carlo  method  to  augment  other  device  and  circuit  simulation  techniques. 

2.2.5  Theory  of  Optical-Phonon  Modes  in  Heterostructures 

During  the  past  year,  we  have  initiated  a  research  project  on  the  theory  of  optical-phonon  modes  in 
heterostructures.  It  has  been  experimentally  verified  that  the  phonon  modes  as  well  as  carrier  wave  func¬ 
tions  are  significantly  <ilfected  by  the  presence  of  heterointerfaces.  Thus,  the  characteristics  of  optical- 
phonon  modes  in  the  layered  structiues  are  distinctively  different  from  their  bulk  counterparts.  Questions 
arise  as  to  the  validity  of  the  conventional  three  dimensional  (bulk)  heatments  of  optical-phonon  model 
in  heterostructure  devices,  especially  when  interactions  by  the  longitudinal-optical  (LO)  phonon  modes 
become  a  dominant  energy  loss  mechanism.  Our  main  objective  is  to  study  the  effects  of  reduced  dimen¬ 
sionality  on  optical-phonon  modes  from  various  aspects. 

In  spite  of  the  short  period  of  time  during  which  this  research  has  been  performed,  we  have  made 
substantial  progress.  Tlircc  papers  have  been  already  published  or  accepted  by  major  journals.  Two  addi¬ 
tional  papers  have  been  submitted  for  publication  and  two  more  are  currently  in  preparation.  As  a  result 
of  the  progress,  this  research  project  will  be  supported  by  separate  grants  from  ONR  and  ARO  which 
have  been  recently  awarded  to  us. 

For  brevity,  this  section  will  focus  on  the  study  of  confined  LO-phonon  frequencies  in  short-period 
strained-layer  superlattices.  Ttie  main  emphasis  will  be  on  the  development  of  a  simple  model  which  can 
calculate  the  frequencies  of  confined  LO  phonons  in  strained-layer  III-V  superlattices  grown  in  the  <001> 
direction.  A  discrete  (rather  than  continuum)  approach  based  on  a  linear-chain  model  is  applied.  The  fre¬ 
quency  shift  arising  simultaneously  from  both  confinement  and  strain  is  accounted  for  by  properly  modi¬ 
fying  the  force  constavxts  to  incorporate  the  strain  effects  into  the  linear-chain  model.  For  a  specific 


Figure  19.  Comparison  of  experimental  response  of  a  GaAs^ 
MSM  photodetector  to  its  simulated  response.  The  finger  width 
is  0.75  Jim  and  the  interfinger  spacing  is  0.5  pm. 


32 


example,  the  frequency  shifts  are  calculated  for  GaAs/GaP  short-period  superlattices,  and  compared  with 
the  existing  Raman  data  of  Armelles  et  al.  [56]  for  the  lowest  order  confined  LO  modes. 


Figure  20  shows  a  typical  schematic  drawing  for  the  linear-chain  model  in  a  superlattice.  As  can  be 
seen  in  the  figure,  the  z  <ixis  is  chosen  as  the  growth  direction  of  the  lattice  in  this  study.  Hence,  for  the 
longitudinal  modes  of  our  interest,  the  atomic  displacement  is  in  the  z  direction.  In  the  absence  of  strain, 
the  force  constant  in  each  layer  is  estimated  based  on  the  frequencies  of  corresponding  bulk  LO  modes  at 
the  zone  center.  Thus,  the  force  constant  in  one  layer  is  modeled  independent  of  the  rest  of  the  structure, 
and  is  the  same  as  the  bulk  value  of  the  material.  The  resulting  characteristic  equations  can  predict  the 
frequency  shifts  due  to  confinement  in  LO  modes  with  accuracy. 

When  a  strain  is  applied,  the  characteristic  equation  of  motion  needs  modification  to  implement  the 
effects  of  strain.  From  the  analyses  of  Ganesan  et  al.  [57]  and  Cerdeira  et  al.  [58],  the  dynamical  equation 
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for  the  LO  modes  with  k=:0  can  be  approximated  as.  to  the  first  order  in  strain. 
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in  bulk  diamond-  and  zinc-blende-typc  semiconductors.  Here,  Uj  is  the  /th  component  of  the  relative  dis¬ 
placement  of  the  two  atoms  in  a  unit  cell:  m  is  the  reduced  mass  of  the  two  atoms;  =  m®o  is  the 
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effective  bulk  force  constant  for  LO  modes  ((a„)  in  the  absence  of  strain;  and  -z —  =  is  the  change 
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in  effective  force  constant  due  to  an  applied  strain  £|,„.  The  indices  i,  k,  1,  and  m  represent  x,  y,  or  z.  In 
Eq.  (39),  the  effect  of  mismatched  lattice  constants  is  included  tlirough  Ein,  which  depends  explicitly  on 
the  lattice  constants.  From  the  symmetry  considerations  in  elasticity  theory,  it  can  be  shown  that  in  cubic 
crystals  the  tensor  K*’’  has  only  tiuce  independent  components  [58]; 
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All  other  components  in  which  x,  y,  or  z  appears  an  odd  number  of  times  become  zero.  Thus,  for  the  LO 
modes  in  the  z  direction,  Eq.  (39)  can  be  simplified  as 

inii,.  =  -  K  <S,’>U2  -  m  {  pe^^  +  q(e,,x  +  Eyy)  }  Uz  -  2mrex2Ux  -  2mr€y2Uy .  (40) 

When  the  strain  is  applied  along  the  major  crystal  axes  (i.e..  x,  y,  and  z),  which  is  the  case  for  strained- 
layer  structures  grown  in  the  <001  >  direction,  the  and  Sy^  terms  vanish  and  the  dynamical  equation 
becomes 
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As  a  result,  the  effect  of  strain  on  the  effective  force  constant  for  LO  modes  can  be  obtained  as 
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Therefore,  within  the  nearest  neighbor  approximation,  the  effects  of  strain  can  be  included  in  the  linear- 
chain  model  for  strained-layer  structures  by  modifying  the  force  constants  in  each  layer  based  on  the  ratio 
given  in  Eq.  (42). 


It  is  well  established  that  the  strain  in  individual  layers  of  a  strained-layer  pseudomorphic  superlat¬ 
tice  can  be  calculated  by  matching  up  lattice  constants  in  the  plane  parallel  to  the  heterointerfaces.  The 
resulting  biaxial  stress  subjects  each  layer  under  biaxial  compressive  or  tensile  strain  in  the  parallel  direc¬ 
tion,  and  uniaxial  tensile  or  compressive  strain  in  the  perpendicular  direction  (i.e.,  along  the  crystal 
growth  direction).  When  isolated  from  the  substrate,  the  lattice  constant  parallel  to  the  interface  an  can  be 
theoretically  obtained  by  minimizing  the  the  elastic  energy  of  each  period  (i.e.,  two  alternating  layers) 
under  the  assumption  that  the  lattice  mismatch  is  accommodated  by  uniform  deformation,  and  is  given  as 
follows  [59]: 
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where  the  index  I  (2)  represents  the  layer  with  smaller  (larger)  bulk  lattice  constant  aj  (a2)  respectively,  d 
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is  the  layer  thickness,  and  f  is  the  mismatch  between  the  two  layers.  The  bulk  shear  modulus  G  can  be 
written  in  terms  of  elastic  constants  C|  |,  Ci2  of  each  material  as  [59]: 
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Then,  the  nonzero  components  of  the  strain  tensor  in  each  layer  are 
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For  the  strained-layer  quantum  wells  and  superlattices  coupled  to  the  substrate,  the  strain  is  accommo¬ 
dated  by  only  one  type  of  layer.  In  this  case,  the  thickness  of  the  other  layer  can  be  considered  very  large 
in  Eq.  (43).  For  comparison  with  other  calculations,  direct  application  of  Eqs.  (42)  and  (45)  to  the  linear- 
chain  model  also  results  in  the  frequency  shift  due  to  strain  in  the  absence  of  confinement  effects  as 
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which  is  correct  to  the  first  order  in  strain,  where  Sj  i  and  S|2  are  elastic  compliances. 

The  frequencies  of  LO-phonon  modes  have  been  calculated  in  the  GaAs/GaP  short-period  strained- 
layer  superlattices  using  the  linear-chain  model  described  above.  The  material  parameters  for  GaAs  and 
GaP  used  in  the  calculation  can  be  found  in  Ref.  60  except  the  p  and  q  values  for  GaAs  which  are  from 
Ref.  61.  Two  different  cases  have  been  considered  for  the  distribution  of  strain:  Table  I  shows  the 
results  when  the  superlattices  are  isolated  from  the  substrates  and,  thus,  the  strain  is  accommodated  by 
both  GaAs  and  GaP  layers  based  on  their  thicknesses  (i.e.,  Eq.  (43)):  and  Table  II  shows  the  results  when 
the  strain  is  applied  only  to  the  G<iP  layers  (i.e..  lattice-matched  to  the  GaAs  substrate).  In  both  cases, 
superlattices  with  two  monolayers  of  GaP  and  two,  four,  and  six  monolayers  of  GaAs  (i.e., 
(GaAs)2/(GaP)2,  (GaAs)^/(GaP).,,  and  (GaAs)^/(GaP)2)  have  been  considered.  The  Am,,  and  Am,  terms 
represent  the  frequency  shifts  due  to  confinement  and  strain  (i.e.,  Eq.  (46))  calculated  in  the  absence  of 
the  other  respectively,  and  the  values  of  Am  are  the  total  shifts  obtained  by  the  model  presented  previ- 


Table  I.  Frequency  shifts  of  confined  LO-phonon  modes  in  GaAs/GaP  short-period  strained-layer  super¬ 
lattices.  The  strain  from  lattice  mismatch  is  accommodated  based  on  the  assumption  that  the  superlattices 
are  isolated  from  the  substrates.  GaAs  stands  for  the  GaAs-like  LO  modes  localized  in  the  GaAs  layers, 
and  GaP  for  the  GaP-like  LO  modes  localized  in  the  GaP  layers.  The  AtOc  and  A©,  terms  represent  the 
shifts  due  to  confinement  and  strain  calculated  in  the  absence  of  the  other,  respectively,  and  the  values  of 
A©  are  the  total  shifts  obtained  by  the  model  presented  in  the  text. 


GaAs  (cm"*) 

A©c+A©, 

A© 

GaP  (cm"*) 
A©c+Aci^ 

A© 

(GaAs)2/(GaP)2 

-4.5 

-5.7 

-22.6 

-21.9 

(GaAs)4/(GaP)2 

1.1 

0.8 

-26.9 

-26.0 

(GaAsy(GaP)2 

1.9 

1.8 

-29.2 

-28.2 

* 


Table  II.  Frequency  shifts  of  confined  LO-phonon  modes  in  GaAs/GaP  short-period  strained-layer 
superlattices.  Ilie  strain  from  lattice  mismatch  is  solely  applied  to  the  GaP  layers.  All  of  the  notations 
are  the  same  as  Table  I.  Therefore,  the  values  of  A©^  for  the  GaAs-like  modes  are  zero  in  this  table. 


GaAs  (cm"*) 

A©c+A©, 

A© 

GaP  (cm"*) 
A©e+A©, 

A© 

(GaAs)2/(GaP)2 

-12.0 

-13.0 

-36.5 

-35.2 

(GaAs)/GaP)2 

-4.1 

4.3 

-36.5 

-35.2 

(GaAs)g/(GaP)2 

-2.0 

-2.1 

-36.5 

-35.3 
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ously. 

It  can  be  seen  from  the  tables  that  the  coupling  effect  ( -  I  Aw  -  (A<»c+A<»s)  I )  is  in  general  small  in 
GaAs/GaP  superlattices.  With  a  3.7  %  lattice  mismatch,  the  largest  deviation  (or  coupling)  is  only  1.3 
cm"'  approximately,  and  decreases  as  the  strength  of  confinement  and  strain  diminishes  in  each  layer  as 
expected.  For  the  GaAs-iike  mode  (i.e.,  the  LO  mode  localized  in  the  GaAs  layer),  six  monolayers  seem 
to  be  thick  enough  to  decouple  the  shifts  due  to  confinement  and  strain  completely.  An  interesting  result, 
however,  is  the  sign  of  deviation.  As  can  be  seen  in  Tables  I  and  II,  the  model  predicts  a  smaller  Aw  for 
the  GaAs-like  mode  (i.e.,  more  negative)  and  a  larger  Aw  for  the  GaP-like  mode  (i.e.,  less  negative)  com¬ 
pared  to  AWc+Aw.;  respectively.  TThis  tendency  can  be  easily  explained  by  a  simple  picture  suggested  by 
the  equivalent  wave-vector  model  [62 1.  Since  the  lattice  constants  are  matched  in  the  plane  parallel  to  the 

interfaces  and  Poisson’s  ratio  has  a  negative  value,  the  lattice  constants  in  the  normal  direction  to  the 

• 

interface  (af*^*  and  af arc  not  equal  and  arc  always  larger  for  the  GaAs  layers  than  for  the  GaP  layers 
(i.e.,  af  >  af Tluis.  the  additive  term  in  the  denominator  of  the  equivalent  wave-vector  should  be 
replaced  by  a  term  smaller  than  1  for  the  GaAs-like  mode  (a°®'Va°“^'’)  and  a  term  larger  than  1  for  the 
GaP-like  mode  (a°“'^Vaf '*'*).  As  a  result,  the  coupling  causes  a  more  negative  Aw  for  the  GaAs-like  mode 
and  a  less  negative  Aw  for  the  GaP-like  mode  compared  to  AWc+AWj,  respectively,  as  mentioned  above. 
Another  interesting  point  is  that  the  confinement  effect  for  the  GaP-like  mode  is  insensitive  to  the  thick¬ 
ness  of  the  GaAs  layers  as  shown  in  Table  II.  Hence,  given  the  distribution  for  strain,  the  phonon  modes 
in  a  pseudomorphic  thin  quantum  well  structure  (not  a  superlattice)  can  be  effectively  studied  by  this 
model.  Unlike  GaA.s/GaP  superlattices.  the  effect  of  coupling  may  be  significant  in  other  material  sys¬ 
tems,  Preliminary  results  show  a  considerably  large  deviation  of  4  cm"'  in  a  (GaAs)2/(InAs)2  superlat¬ 
tice. 

The  calculated  data  for  frequency  shift  Aw  have  been  compared  with  the  Raman  measurements  of 
Armelles  et  al.  [56].  Tlie  results  for  the  GaP-like  mode  and  the  GaAs-like  mode  are  plotted  in  Fig.  21  and 
Fig.  22  respectively.  In  these  figures,  LO„  represents  the  ;jth  order  confined  LO  mode.  The  data  for  the 
(GaAs)^/(GaP).,  superlattice  in  Tables  I  and  II  are  denoted  as  #1  of  (a)  and  #2  of  (a)  respectively.  In 
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addition,  the  results  for  a  (GaAs)^/(GaP)^  superlattice  are  included  with  two  different  configurations  for 
the  accommodation  of  strain:  #1  of  (b)  shows  the  results  when  the  superlattice  is  isolated  from  the  sub¬ 
strate  (which  corresponds  to  #1  of  (a)):  and  #2  of  (b)  shows  the  results  with  a  0.7  %  strain  in  the  GaAs 
layers  and  a  3.0  %  strain  in  the  GaP  layers.  The  strain  distributions  used  to  obtain  data  #2  of  (a)  and  #2  of 
(b)  are  identical  to  those  used  in  the  <uialysis  of  Armelles  et  al.;  that  is,  for  #2  of  (a)  all  of  the  3.7  %  strain 
is  in  the  GaP  layers  and  for  #2  of  (b)  the  strain  distribution  is  as  described  previously.  As  can  be  seen 
from  the  figures,  the  calculated  data  arc  in  good  agreement  with  the  Raman  measurements.  Indeed,  the 
peaks  in  the  Raman  spectra  are  well  placed  between  the  calculated  results  for  the  two  different  strain 
configurations.  This  level  of  agreement  suggests  that  the  measured  frequency  shift  can  be  explained  by 
this  linear-chain  model  with  a  realistic  configuration  of  strain  accommodation.  Based  on  the  overall  com¬ 
parison  between  theoretical  predictions  and  experimental  data,  it  is  concluded  that  the  only  available 
Raman  frequency  shift  data  [56]  on  GaAs/GaP  short-period  strained-layer  superlattices  are  better 
explained  by  assuming  that  the  superlattices  are  isolated  from  the  substrates  than  by  assuming  that  they 
are  lattice  matched  to  the  substrates.  While  this  conclusion  has  not  been  reached  in  previous  analyses,  it 
is  not  imreasonable  considering  the  critical  thickness  for  the  formation  of  misfit  dislocation  near  the 
substrate-superlattice  interface.  An  alternative  explanation  is  that  the  Raman  measurements  by  Armelles 
et  al.  [56]  characterized  regions  of  the  GaAs/GaP  superlattices  that  were  only  partially  decoupled  from 
the  substrate.  However,  for  a  more  detailed  comparison  between  experimental  and  theoretical  data, 
further  clarification  is  essential  concerning  experimental  uncertainties  such  as  the  depth  of  Raman  probe, 
inhomogeneities  in  the  strain  distribution  along  the  growth  direction,  and  interface  roughness. 

An  effort  has  been  made  to  investigate  higher  modes  as  well.  For  the  (GaAs)g/(GaP)2  superlattice, 
the  theoretical  calculation  for  the  GaAs-like  LOi  modes  predicts  -  287  cm"'  (when  isolated  from  the 
substrate)  and  ~  283  cm"'  (when  lattice  matched  to  the  substrate),  respectively.  The  experimental  result 
is  -  290  cm"'  which  shows  good  agreement  with  the  theoretical  calculations.  The  Raman  data  for  the 
GaP-like  higher  modes  in  the  (GaAs)^/(GaP).j  superlattices  were  not  available  for  comparison.  Due  to 
the  simplicity  of  our  model,  it  is  expected  that  the  theoretical  prediction  for  the  higher  modes  deteriorates 
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rapidly  as  the  corresponding  equivalent  wave-vector  increases. 

2.2.6  Simulation  of  Delta-Doped  FETs 

This  work  has  focused  on  electron  transport  and  device  performance  in  delta-doped  field  effect 
transistor  (FET)  structures.  The  superior  conduction  property  of  delta-doped  FETs  over  that  of  conven¬ 
tional  (i.e.,  uniformly  doped)  FET  structures  makes  this  device  one  of  the  most  promising  candidates  for 
ultra-high  speed  digital  and  millimeter-wave  applications.  High  electron  mobility,  reduced  DX  center 
concentration,  improved  threshold  controllability,  as  well  as  high  breakdown  voltage  are  important 
advantages  of  delta-doped  FETs.  Studies  of  electron  transport  and  electrical  characteristics  for  such  dev¬ 
ices  are  both  necessary  and  important. 

Our  approach  to  this  research  project  is  two-fold;  we  employ  ensemble  Monte  Carlo  simulations 
for  submicron  AlGaAs/GaAs  delta-doped  HEMTs  and  two-dimensional  drift-diffusion  simulations  foP 
delta-doped  GaAs  MESFETs.  The  objective  of  Monte  Carlo  simulations  is  to  provide  physical  insight 
into  the  detailed  operation  of  submicron  delta-doped  HEMTs.  Scattering  properties,  real  space  transfer 
and  field  dependence  of  active  device  dimensions,  doping  concentration  and  bias  conditions  are  studied  to 
provide  better  understanding  of  device  operation.  A  two-dimensional  drift-diffusion  model  is  used  to 
analyze  the  electrical  characteristics  of  a  delta-doped  GaAs  MESFET,  providing  a  general  guide  for  more 
practical  device  design.  Currently,  we  have  finished  the  program-development  stage  and  have  started  to 
obtain  and  analyze  preliminaiy  data.  Two  of  papers  will  be  submitted  to  the  journals  in  the  near  future. 
We  also  plan  to  submit  a  paper  to  the  1990  lEDM  which  will  be  held  in  December.  The  models  used  in 
the  simulations  and  the  preliminary  data  vvill  be  discussed  briefly  as  follows: 

A.  Ensemble  Monte  Carlo  Simulation 

The  ensemble  Monte  Carlo  program  is  developed  to  model  transport  properties  of  electrons  in 
heterojunction  FET  stnictures.  In  this  program,  which  incorporates  a  non-parabolic  F-L-X  analytical 
band  structure,  the  following  scattering  mechanisms  are  included:  polar  optical  phonon,  iatervalley,  ion- 
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ized  impurity,  electron-electron,  and  impact  ionization.  Electrostatic  potential  and  the  electric  field  distri¬ 
bution  within  a  device  are  evaluated  for  each  field-adjusted  time  interval  by  solving  the  two-dimensional 
Poisson  equation  self-consistently.  The  effect  of  carrier  degeneracy  is  considered  by  employing  Fermi 
statistics  with  calculated  local  quasi-Fermi  level  and  energy  distribution.  Real  space  transfer  of  electrons 
between  different  materials  is  also  included  in  the  program.  A  weighted  injection  scheme  is  used  to 
improve  simulation  efficiency.  Also,  an  attempt  is  made  to  take  the  effect  of  surface-defect-related 
Fermi-level  pinning  into  account.  This  program  is  capable  of  simulating  both  delta-  and  uniformly-doped 
HEMT  structure  (hereafter,  D-HEMT  and  U-HEMT,  respectively).  Theoretical  calculations  and  com¬ 
parisons  are  performed  for  AlGaA.s/GaAs  D-HEMTs  and  U-HEMTs.  The  superior  performance  of  the 
D-HEMT  in  terms  of  device  structure,  doping  and  bias  condition  is  demonstrated.  Comparison  with 
available  experimental  results  will  be  made  to  verify  this  theoretical  study. 

Distribution  of  electron  concentration,  field,  drift  velocity,  energy,  and  valley  occupancy  have  beefc 
calculated  for  AlGaAs/GaAs  D-HEMTs  and  U-HEMTs  with  different  device  featine  sizes,  doping,  and 
bias  conditions.  Comparison  of  electrical  characteristics  between  the  two  device  structures  is  made  based 
on  these  calculations.  In  order  to  make  the  comparison  on  a  one-to-one  basis,  devices  with  the  same  chan¬ 
nel  length  and  under  the  same  bias  are  simulated.  Doping  of  both  devices  are  such  that  they  have 
equivalent  total  charge  when  the  impurities  are  fully  ionized.  Two  specified  conditions  are  chosen  for 
study: 

(1)  both  devices  have  same  distance  from  gate  to  channel: 

(2)  both  devices  have  same  tlueshokl  voltage. 

The  parameters  used  in  the  simulatiotis  are  summarized  in  Table  III 

>From  our  preliminary  results,  we  can  draw  the  following  conclusions: 

•  For  both  conditions  (1)  and  (2).  the  channel  electron  density  of  D-HEMT  is  higher  than  that  of  U- 
HEMT.  For  U-HEMT  in  condition  (1),  this  is  due  to  the  fact  that  comparatively  higher  field  is  found  in 
the  channel  near  the  drain  edge  of  the  gate.  High  electric  field  causes  local  electron  heating,  and  heated 
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electrons  with  sufficient  energy  will  experience  intcrvalley  scattering  and  real  space  transfer.  Simulation 
results  show  that  there  are  more  real  space  transfer  events  and  intervalley  scatterings  in  this  case.  For  U- 
HEMT  in  condition  (2).  decreased  electron  density  in  the  channel  is  caused  by  the  low  doping  in  the 
AlGaAs  layer  and  the  longer  distance  from  gate  to  channel. 

•  For  both  conditions.  Iiigher  current  and  higher  transconductance  are  obtained  for  the  D-HEMT. 
This  is  partly  because  of  the  larger  channel  electron  density  in  D-HEMTs,  as  discussed  above.  In  addi¬ 
tion,  higher  drift  velocity  as  a  result  of  relatively  smaller  electric  field  in  the  D-HEMT  channel  (fewer 
velocity  randomizing  collisions)  further  improves  device  performance.  This  trend  is  in  accord  with  a 
macroscopic  consideration  from  a  simple  analytical  model.  Analytical  and  measurement  results  show 
that  the  threshold  voltage  for  D-HEMT  and  U-HEMT  can  be  expressed  as: 

D-HEMT:  V„,  = -qNjdfw  -  0.5d)/8  (47) 

» 

U-HEMT :  V,,,  =  -qNdW-/(2£) 

Thus,  for  parameters  given  in  condition  (1).  high  V,i,  is  expected  for  the  U-HEMT,  which  in  turn  results 
in  smaller  current  and  transconductance  for  a  given  bias  condition.  For  condition  (2),  longer  distance  d 
from  the  gate  to  the  channel  should  be  used  in  order  to  have  the  same  threshold  voltage.  However,  this 
deceases  the  controllability  of  gate  voltage  on  channel  electrons  and  causes  decreased  transconductance. 
Within  the  Gradual  Channel  Approximation  (GCA),  the  saturation  current  and  transconductance  are 
inversely  proportional  to  the  gate-channel  distance.  Thus,  increasing  this  distance  causes  degraded  per¬ 
formance. 

•  Simulation  results  show  that  as  gate  voltage  increases  above  0  V  (both  devices  are  in  depletion 
mode),  the  transconductance  decreases  even  though  the  channel  field  is  small.  This  suggests  parallel  con¬ 
duction  paths  exist  in  the  doped  AlGaAs  layer,  where  field  is  much  higher  and  much  more  impurity 
scatterings  are  expected.  This  effect  is  more  severe  in  condition  (1)  for  U-HEMT  since  the  electric  field  in 
the  AlGaAs  layer  increases  greatly  as  doping  increases  and  distance  decreases. 
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Our  results  demonstrate  the  superior  conducting  property  of  D-HEMTs  over  U-HEMTs.  Another 
feature  which  further  degrades  U-HEMTs  with  doping  over  1.0x10*®  cm“®  is  due  to  the  presence  of  a 
deep-level  trap  (DX  center)  in  the  AlGaAs  layer,  which  is  believed  to  be  less  prevalent  in  the  D-HEMT 
structures. 

B.  Two-Dimensional  Drift-Diffusion  Sinntlarion 

In  this  program,  coupled  phenomenological  semiconductor  equations  are  solved  numerically  to 
simulate  electrical  characteristics  of  delta-doped  and  conventional  GaAs  MESFETs.  These  equations  are: 
Poisson  equation,  current  continuity  equation,  and  current  density  equation.  A  mcdified  two  piece 
velocity-field  relation  is  employed  to  lake  current  saturation  into  account.  The  Fermi  distribution  is  used 
in  the  calculation.  A  Gaussian  doping  profile  is  generated  to  fit  measurement  of  experimental  profiles 
resulting  from  ion  implantation.  Distributions  of  electron  concentration,  potential,  field,  and  current  de% 
sity  are  simulated  for  various  operating  conditions.  Small  signal  parameters  (capacitance,  transconduc¬ 
tance)  are  determined  to  predict  device  performance.  Comparison  is  made  of  delta-doped  and  conven¬ 
tional  uniformly-doped  MESFET  (hereafter,  D-MESFET  and  U-MESFET,  respectively)  calculations  with 
experimental  results.  Using  a  Crank-Nicolson  scheme,  the  program  is  able  to  simulate  transient  perfor¬ 
mance.  Based  on  this  model,  current-voltage  and  capacitance-voltage  characteristics  can  be  easily  calcu¬ 
lated.  Good  agreement  between  simulation  and  measurement  is  obtained.  Parameters  used  in  calculation 
and  from  mcasurcmental  results  arc  summarized  in  Tabic  IV. 

Our  preliminary  results  suggest  the  following  conclusions: 

•  Channel  length  L  and  gate-to-channel  distance  d  arc  the  two  most  important  dimensional  parame¬ 
ters.  The  tran.sconductancc  increases  as  L  and  d  decrease.  However,  a  decrease  in  gate-to-chaimel  dis¬ 
tance  increases  the  channel  field  and.  thus,  deteriorates  the  breakdown  voltage.  Spacer  distance  from  gate 
to  drain  is  crucial  in  determining  the  maximum  electric  field,  especially  for  self-aligned  structures.  By 
increasing  this  distance,  the  maximum  iicid  can  be  reduced.  However,  at  the  same  time  the  parasitic 
series  resistance  increases  which,  in  turn,  reduces  the  transconductance.  Half  maximum  doping  width  is 


also  an  important  parameter  in  terms  of  transconductance.  Trade-off  should  be  made  among  these  param¬ 
eters  in  device  design  in  order  to  get  optimized  performance. 

•  Increasing  maximum  (or  peak)  doping,  in  D-MESFET  increases  electron  concentration.  This 
is  the  main  contribution  to  the  high  current  capability  of  the  structiue.  However,  as  Nn,,^  increases,  the 
maximum  field  increases  as  well.  This  will  cause  reduced  drift  velocity  and  sacrifices  the  transconduc¬ 
tance  as  reaches  to  a  critical  point.  A  slight  increase  of  background  doping  in  the  cap  layer  is  useful 
in  reducing  the  electric  field  in  the  channel.  For  D-MESFETs,  where  the  conducting  channel  is  a 
heavily-doped  layer,  a  decreased  field  in  this  layer  improves  saturation  performance. 

•  Compared  with  conventional  the  U-MESFET,  the  D-MESFET  shows  improved  performance  in 
cunent  capacity  and  transconductance.  The  most  important  feature  of  D-MESFET  is  its  very  high  chan¬ 
nel  carrier  concentration  (factor  of  10  higher  than  U-MESFET).  Generally,  electrons  in  such  a  heavily- 
doped  region  undergo  more  impurity  .scatterings  and  exhibit  low  mobility.  However,  at  high  temperature 
and  high  electric-field,  from  which  the  electron  can  gain  energy,  impurity  scattering  becomes  less 
significant  on  overall  mobility.  Moreover,  when  the  electric  field  exceeds  a  critical  value,  the  saturated 
velocity  determines  device  performance.  Both  theoiy  and  experimental  results  show  that  saturated  drift 
velocity  does  not  depend  on  the  low  field  mobility,  i.e.  the  same  saturation  velocity  for  D-MESFET  and 
U-MESFET  is  expected.  Another  feature  which  makes  the  D-MESFET  performance  superior  to  that  of 
the  conventional  U-MESFET  is  that  for  the  same  gate-to-channel  distance,  lower  maximum  electric  field 
is  expected  because  of  the  linear  conduction  band  edge  dependence  on  gate-to-channel  distance.  In  a  con¬ 
ventional  U-MESFET,  a  quadratic  dependence  of  the  conduction  band  edge  on  the  gate-to-channel  dis¬ 
tance  is  expected  due  to  its  constant  doping  concentration  under  the  gate.  Reduced  field  improves  satura¬ 
tion  and  breakdown  pcrfonnance. 

2.3  Publications  and  Presentations 

During  the  last  year  we  have  made  seven  oral  presentations  at  national  and  international  confer¬ 
ences.  Also,  twelve  written  manuscript.^  have  been  published  in  the  refereed  literature.  In  addition,  three 
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Other  manuscripts  have  been  accepted  for  publication  and  eight  manuscripts  have  been  submitted  or  are 
currently  in  preparation.  The  following  paragraphs  summarize  the  presentations  and  publications  made 
under  this  program  during  the  last  year,  and  include  material  not  previously  reported  to  ONR  in  the  last 
renewal  proposal  as  well  as  material  which  was  accepted  for  publication  in  1989  and  has  since  been  pub¬ 
lished  in  1990. 
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A.  Conference  Presentations  and  Seminars 

D,L.  Woolard,  J.L.  Pelouard,  R.J.  Trew,  M.A.  Littlejohn  and  C.T.  Kelly,  "Hydrodynamic  Hot-Electron 
Transport  Simulations  Based  on  The  Monte  Carlo  Method",  presented  at  the  Sixth  International  Confer¬ 
ence  on  Hot  Carriers  in  Semiconductors,  July  23-28, 1989,  Scotttsdale,  AZ. 

L.F.  Register,  M.A.  Littlejohn  and  M.A.  Stroscio,  "Polar  Optical  Phonon  Scattering  of  Charge  Carriers  in 
Alloy  Semiconductors;  Effects  of  Phonon  Localization",  presented  at  the  Sixth  International  Conference 
on  Hot  Carriers  in  Semiconductors,  July  23-28, 1989,  Scottsdale,  AZ, 

J. L,  Pelouard.  M.A.  Littlejohn  and  H.P.  Belgal,  "Monte  Carlo  Study  of  Ballistic  and  Quasi-Ballistic  Elec¬ 
tron  Transport  in  Semiconductors,"  presented  at  the  Sixth  International  Conference  on  the  Numerical 
Analysis  of  Semiconductor  Devices  and  Integrated  Circuits,  July  11-14, 1989,  Dublin,  Ireland. 

L, F.  Register.  M.A.  Littlejohn  and  M.A.  Stroscio.  "Path-Integral  Monte  Carlo  Methods  for  Ultra-Small 
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as  the  Associate  Dean  of  Academic  Affairs  in  the  College  of  Engineering.  Dr.  Hauser  has  recently 
assumed  major  technical  and  administrative  responsibilities  with  the  newly  -  formed  NSF  Engineering 
Research  Center  for  Advanced  Electronic  Materials  Processing.  Both  these  former  investigators  on  this 
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Drs.  Littlejohn  and  Kim  ser\'e  as  co-principal  investigators  on  this  research  effort.  Their  resumes 
are  in  Appendix  C.  Tliey  are  responsible  for  day-to-day  management  and  direction  of  the  research  pro¬ 
gram  aspects.  Three  graduate  students  have  been  supported  on  this  program,  and  we  plan  to  add  a  fourth 
student  during  the  fall  semester. 
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Simulation  results  of  the  temporal  evolution  of  photocurrent  in  an  interdigitated  GaAs  metal- 
semiconductor-meta!  photodetector  are  presented.  The  dependence  of  response  time  on  the 
distance  between  fingers  (0. 1  and  0.25  ptm)  is  investigated.  The  solutions  of  the  time- 
dependent  Schrodinger  equation  and  ensemble  Monte  Carlo  calculations  are  employed.  For  a 
device  with  0. 1  ptxa  finger  spacing,  the  response  time  of  an  intrinsic  detector  is  less  than  2  ps, 
with  hole  current  decay  being  the  major  limiting  factor.  The  role  of  parasitics  is  shown  to 
significantly  increase  the  simulated  response  time. 


It  has  recently  been  shown  that  GaAs  metal-semicon- 
ductor-metal  (MSM)  photodetectors' ’  can  be  used  as  high 
performance  optical  detectors  in  the  0.8-0.87  /im‘ wave¬ 
length  range.  High  bandwidth  monolithic  optoelectronic  re¬ 
ceivers  using  MSM  photodetectors  have  already  been  dem¬ 
onstrated.*"'®  Because  of  the  demand  for  even  faster  circuits, 
there  is  a  need  for  numerically  simulated  results,  especially 
for  detectors  with  dimensions  even  smaller  than  those  con¬ 
sidered  previously.’ 

In  this  letter  we  report  theoretical  results  on  the  dynam¬ 
ic  behavior  of  photogenerated  carriers  in  GaAs  MSM  photo¬ 
detectors  with  sub-half-micron  finger  spacings.  The  time 
evolution  of  photocurrent  is  calculated  for  detectors  with  a 
0. 1  and  0.25  /rm  distance  between  metal  fingers.  The  detri¬ 
mental  influence  of  the  extrinsic  circuit  elements  on  the  de¬ 
vice  speed  is  also  shown. 

The  GaAs  MSM  photodetector  is  shown  schematically 
in  Fig.  1(a).  The  device  is  modeled  as  a  one-dimensional 
structure  having  an  active  layer  of  length  L  with  two  ideal 
Schottky  barrier  contacts.  The  background  doping  level  in 
the  active  layer  was  chosen  to  be  5x  lO''*  cm"-’  (n  type). 
Direct-write  electron  beam  lithography  can  be  used  to  ex¬ 
pose  such  a  sub-half-micron  interdigitated  finger  pattern. 
Because  of  the  low  doping  concentration,  the  optically  sensi¬ 
tive  regions  between  the  fingers  are  completely  depleted  even 
at  very  small  bias  voltages.  Thus  we  may  draw  the  energy- 
band  diagram  as  shown  in  Fig.  1(b).  Under  somewhat  sim¬ 
plifying  conditions,  we  assume  that  carriers  are  generated 
halfway  between  the  fingers  with  a  Gaussian  spatial  distribu¬ 
tion. 

We  have  employed  two  different  methods  to  predict 
time  responses.  A  quantum  mechanical  analysis  uses  nu¬ 
merical  solutions  of  the  time-dependent  Schrodinger  equa¬ 
tion  in  the  effective  mass  approximation.  The  self-consistent 
Monte  Carlo  simulation  us^  in  this  study  is  described  else¬ 
where,"  and  employs  a  set  of  GaAs  parameters  published 
previously.'*  Neither  quantum  effects  (tunneling)  nor  ex- 
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ternal  parasitic  effects  (capacitance)  have  been  included  in 
this  simulation.  The  boundary  conditions  assume  the  conti¬ 
nuity  of  the  quasi-Fermi  level  through  the  metal-semicon¬ 
ductor  interface  for  each  carrier  under  reverse  bias.  Because 
the  barrier  heights  for  electrons  and  holes  and  re¬ 
spectively)  are  much  larger  than  kT,  the  current  densities 
due  to  injected  carriers  from  the  metal  into  the  semiconduc¬ 
tor  are  negligible. 

The  solutions  of  the  time-dependent  Schrodinger  wave 
equation  for  an  electron  and  a  hole  are  shown  in  Fi^.  2(a) 


MSM 


FIG.  1.  (a)  Schematic  view  of  an  interdigitated  GaAs  MSM  photodetector. 
The  distance  between  lingers  is  equal  to  £.  (b)  Energy  band  diagram  of  a 
detector  under  bias  i  The  initial  ( r  ==  0)  distnbuttons  of  photocarriers  are 
Gaussian  centered  i  the  middle  of  the  device.  The  mitial  energies  of  elec¬ 
trons  and  holes  arc  -  *  and  1  meV,  respectively. 
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Analysis  of  a  GaAs  Metal-Semiconductor-Metal 
(MSM)  Photodetector  with  0. 1-fim 
Finger  Spacing 

WACLAW  C.  KOSCIELNIAK.  MICHAEL  A.  LITTLeOHN,  member,  ieee,  and  JEAN-LUC  PELOUARD 


Abstract— Tht  temporal  response  of  a  GaAs  raetal-semiconductor- 
metal  (MSM)  photodelector  with  a  finger  spacing  of  0.1  fim  has  been 
analyzed.  The  intrinsic  detector  has  been  found  to  have  a  minor  effect 
(25-pcrccnt  increase)  on  the  full  width  at  haif^maximum  (FWHM)  of  the 
temporal  response  of  the  device  and  its  parasitic  circuit  elements.  The 
analysis  indkates  that  a  long  lime  constant  due  to  the  decay  of  holes  is 
solely  responsiUe  for  this  increase.  The  smallest  FWHM  for  this  detector 
is  estimate  to  be  less  than  2.5  ps. 

I.  Introduction 

Applications  of  a  GaAs  metal-semiconductor-metal 
(MSM)  photodetector  have  been  demonstrated  by  several 
researchers  (i]-(5].  Recently,  photodetectors  with  a  very  short 
response  time  and  a  high  quantum  efficiency  have  been 
realized  (3],  [6).  The  0.8-0.87-#tm  near-infrared  wavelength 
range  of  GaAs  detectors  has  been  extended  to  1.3-l.S  fim  in 
InGaAs-based  devices  [7],  [8].  In  addition,  other  ways  of 
increasing  the  detector's  responsivity  have  been  investigated 
[9].  The  availability  of  high-performance  optical  detectors 
compatible  with  MESFET  processing  make  the  MSM  device 
an  ideal  choice  for  monolithically  integrated  optoelectronic 
circuits  [10]-(14]. 

To  our  knowledge,  there  are  few  reports  on  simulation 
results  for  this  device  (61,  (15].  In  this  letter,  we  present  a  one¬ 
dimensional  analysis  of  the  temporal  response  of  a  GaAs  MSM 
photodetector  with  sub-haif-micrometer  finger  spacing  (0.1 
/tm).  An  intrinsic  device  with  low-level  optical  excitation  and 
different  electric  field  profiles  is  studied.  We  also  include  an 
extrinsic  circuit  to  determine  the  importance  of  parasitic 
circuit  elements  on  the  response  speed  of  this  device. 

II.  Analysis  and  Results 

The  photodetector  structure  is  schematically  shown  in  Fig. 
1.  The  unintentionally  doped  epitaxial  layer  (n-lype)  has  a 
background  donor  concentration  of  5  x  10'^  cm"’.  The  use  of 
tungsten  silicide  (WSi^)  for  Schottky  contacts  would  make  this 
device  fully  compatible  with  a  standard  refractory-gate  MES- 
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Fig.  1 .  Structure  of  an  interdigitated  GaAs  MSM  phoiodetector.  The  device 
consists  of  30  0.l-/jm  fingers  and  the  photosensitive  area  is  5.9  x  8  /im^. 

FET  process  (13).  Direct-write  electron-beam  lithography  can 
be  used  to  expo-  'he  multifinger  pattern. 

We  use  a  one  .nensional  model  to  obtain  estimates  of  the 
minimum  time  ponse  of  this  device.  In  order  to  simulate 
potential  effect  -ithin  the  intrinsic  detector,  we  introduce 
spatial  variation  •  f  acceptor  or  donor  concentrations  to  create 
two  built-in  electric  field  profiles,  which  differ  from  those  of 
the  undoped  device.  The  doping  profiles  can  be  created,  for 
example,  by  focused-ion-beam  implantation  of  the  interfinger 
regions,  as  illustrated  in  Fig.  2(a).  From  Poisson’s  equation, 
we  find  the  electrostatic  potential  variations  and  the  resulting 
energy-band  diagram.  For  no  illumination,  the  conduction 
band  edges  for  the  three  cases  considered  here  are  plotted  in 
Fig.  2(b)  for  an  applied  voltage  of  0. 1  V  between  adjacent 
fingers.  Due  to  the  very  low  doping  concentration  of  the  active 
layer,  our  calculations  reveal  that  the  regions  between  fingers 
are  completely  depleted  of  free  carriers  at  0. 1  V  for  all  three 
doping  profiles,  and  this  would  remain  true  if  the  finger 
pattern  were  processed  directly  on  a  semi-insulating  substrate. 
The  electric  field  is  constant  without  dopant  modification  (see 
the  middle  curve  in  Fig.  2(b)).  With  the  acceptor-doping 
profile  (upper  curve),  an  increased  electric  field  occurs  at  the 
positively  biased  electrode  (right-hand  boundary)  forcing  the 
photoexcited  electrons  to  move  faster  in  that  region,  relative  to 
the  undoped  case.  Conversely,  with  the  donor-doping  profile 
(bottom  curve),  a  higher  electric  field  close  to  the  negatively 
biased  electrode  forces  the  photoexcited  holes  to  move  faster 
in  that  region,  relative  to  the  undoped  case.  Intuitively,  it 
would  appear  that  the  last  of  these  three  cases  could  be  the 
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The  ability  to  model  equilibrium  electronic  properties  of  ultrasmall  devices  us* 
ing  a  path-integral  Monte  Carlo  (PIMC)  method  is  demonstrated.  First,  a  direct 
sampling  PIMC  method  and  its  advantages  over  Metropolis  importance  sampling 
PIMC  methods  in  this  application  are  described.  Then  two  potential  structures 
typical  of  ultrasmall  devices,  a  single  and  a  coupled  double  finite  square  well  poten¬ 
tial,  are  analyzed  and  the  results  compared  to  exact  analytic  results.  The  latter 
example  demonstrates  the  ability  to  accurately  resolve  coupling  induced  energy 
level  splittings  at  least  as  small  as  2%  of  the  nominal  energy.  (Lower  levels  of 
coupling  were  not  considered.)  An  additional  example  given,  for  which  analytic 
results  are  not  available,  is  the  analysis  of  a  two-dimensional  nfinite,  periodic 
quantum  well  array.  The  PIMC  technique  described  here  mak  -  way  for  the  ap¬ 
plication  of  the  Feynman  path-integral  formalism,  an  advantax*-  us  formalism  for 
treating  scattering  processes,  to  the  study  of  the  electronic  propr'iies  of  ultrasmall 
devices. 


1  Introduction 

Proper  treatment  of  scattering  processes  such  as 
I  charge-carrier-phonon  coupling  is  crucial  for  accurate 
j  modeling  of  the  electronic  properties  of  semiconduc- 
^  tor  devices.  A  powerful  and  conceptually  appealing 
treatment  of  carrier-phonon  coupling  is  provided  by 
the  Feynman  Path-Integral  (FPI)  formalism  of  quan¬ 
tum  mechanics.*’*  The  power  of  the  FPI  formalism 
in  treating  coupling  processes  was  demonstrated  6rst 
in  quantum  electrodynamics  (QED),*  but  soon  its  ad¬ 
vantages  in  the  formally  similar  treatment  of  carrier- 
phonon  coupling  in  semiconductors  also  wu  realized^ 
and,  since,  several  analytic  FPI-based  studies  of  the 
electronic  properties  of  bulk  Kmiconductors  have  fol¬ 
lowed  (for  example,  References  5-10  and  references 


therein).  In  recent  years,  advances  in  semiconduc¬ 
tor  technologies,  have  led  to  the  realisation  of  devices 
whose  spatid  dimensions  are  on  the  same  Kale  as  the 
thermal  deBroglie  wavelength  of  the  charge  carriers. 
TheK  devices  exhibit  behavior  that  is  either  influenced 
or  governed  by  quantum  interference  effects,  and  they 
must  be  modeled  accordingly.  Further,  on  the  ultra- 
small  spatial  and  temporal  scales  at  which  these  devices 
operate,  the  Kmiclassical  treatment  of  carrier-phonon 
scattering  u  spatially  localised,  instantaneous  events  is 
of  reduced  value.  However,  despite  the  fully  quantum 
mechanical  nature  of  the  FPI  formalism,  there  has  not 
been  widespread  application  of  it  to  the  study  of  the 
electronic  properties  of  these  ultrasmall  structures,  De¬ 
spite  it's  conceptual  advantages,  application  of  the  FPI 
formalism  poses  significant  pr,  .tical  problems.  For  an- 
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ABSTRACT 


The  Heterojunction  Bipolar  Transistor  (HBT)  is  not  a  new  device.  In  the  late  1940's  W.  Shockley^ 
gave  the  basic  ideas  and,  about  ten  years  later,  H.  Kroemer^  published  the  first  paper  to  show  the 
praising  advantages  of  this  new  structure.  The  first  device  applications  have  been  demonstrated  in 
AlGaAs/GaAs  materials,  taking  advantage  of  the  natural  lattice  match  in  this  system.  Currently,  high* 
performance  G^s-based  HBTs  have  been  reported  and  the  InP-based  transistors  are  still  promising.  In 
this  paper  we  will  review  the  reasons  why  this  device  has  required  more  than  forty  years  of  development 
and  what  its  future  can  be  in  the  InP-based  materials  systems. 


1.  INTRODIJCTTON 

In  bipolar  transistors  the  main  current  is  carried  by  electrons  injected  into  the  base  by  the  emitter-base 
junction  and  collected  by  the  base-collector  junction.  There  is  a  transistor  effect  if  the  collector  current  is 
nearly  equal  to  the  emitter  current.  Unfortunately,  electron  injection  (J„)  into  the  base  accompanies  hole 
injection  (jp  into  the  emitter  creating  an  additional  parasitic  current.  Because  the  current  densities  in  the 
homojunction,  or  Classical  Bipolar  Transistor  (CBT),  are  given  by  the  Shockley  formula  for  short  base 
diodes,  the  injection  ratio  (Jp/J„)  depends  on  the  ratio  of  base  (N^)  to  emitter  (Nq)  doping  levels: 


(1) 


:n  order  to  obtain  an  injection  ratio  small  enough  for  significant  current  gain,  the  CBT  must  have  a  much 
higher  doping  level  in  the  emitter  than  in  the  base.  Thus,  the  CBT  will  have  a  high  resistivity  base.  This 
has  various  well  known  negative  consequences  on  the  behavior  of  the  transistor  (Early  effect,  reduction  of 
frnax*  etc.). 


In  contrast,  in  the  first  approximation,  the  energy  discontinuity  in  the  HBT  emitter-base  heterojunction 
introduces  an  additional  term  to  the  injection  ratio  equation: 


(2) 
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Monte  Carh^generated  velocity  diitrihuiion  functions  have  been  utiliied  to  develop  criteria  to 
elauify  electrons  eonirihuting  to  non* stationary  transport  in  semiconductors.  These  classifica* 
tions  are  given  in  terms  of  ballistic,  quasi* ballistic  and  velocity*rtlaxed  electron  populations.  An 
analytical  model  has  been  built  which  describes  these  populatioiu.  While  the  results  Aave  been 
0 6 tamed  for  GaAs  planar* doped  barriers,  we  believe  that  the  definitions  and  concepts  «ne  va!id 
in  a  more  general  sense.  To  support  this  proposition,  some  aspects  of  lnP*bastd  heterojunetion 
6ipo/ar  tranjiitori  will  also  be  discussed. 


1  Introduction 

In  1985,  two  group!  reported  experimental  evidence  (or  hot  electron  transport  in  Larrier 
injection  devices  (l,2l.  These  and  later  results  established  the  existence  of  a  peak  in  the  elec¬ 
tron  energy  distribution  function  at  the  barrier  injection  energy  [3,4].  This  peak,  observable  a 
significant  distance  from  the  barrier,  was  attributed  to  ballistic  transport  in  both  homojunction 
planar-doped  barriers  and  heterojunetion  barriers* 

Recently,  Hess  and  lafrate  explored  these  results  in  an  effort  to  elucidate  further  the  pro¬ 
cesses  of  ballistic  transport  and  velocity  overshoot  (5).  They  developed  definitions  for  these  terms, 
clusified  the  basic  experimental  results,  and  concluded  that,  while  the  basic  physics  of  ballistic 
transport  is  well  understood  qualitatively,  much  quantitative  work  remains  to  be  done,  especially 
with  respect  to  device  applications. 

The  purpose  of  this  paper  is  to  investigate  transport  in  the  base  region  of  GaAs  planar- 

doped  barriers  used  as  ballistic  launchers.  Our  results  expand  the  definitions  of  Hess  and  lafrate 

and  further  classify  electrons  in  the  velocity  distribution  function  as  to  their  role  In  ballistic  and 

quasi-ballistic  transport  in  semiconductors. 
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Tht  application  of  Path* Integral  Monte  Carlo  (PIMCJ  methodi  to  the  study  of  the  electronic  prop* 
erties  of  ultrasmall  devices  is  du cussed.  First,  as  an  introduction,  the  advantages  of  the  Feynman 
path* integral  (FPJ)  formalism  for  this  app/tcation  are  desenfred.  Nett  the  FPI  formalism  itself 
is  reviewed  briefly ,  Then^  a  practical  PIMC  method  for  modeling  equilibrium  conditions  in  these 
devices  IS  iesert^ed.  last,  lacking  a  nonequilibrium  PIMC  metAod  proven  widely  applicable  to 
ulirasmall  device  modeling ^  vanous  basic  approaches  and  their  /imitations  are  descn'Aed  to  ethibit 
the  fundamental  problems  yet  to  be  overcome. 


1  Introduction 

Proper  treatment  of  scattering  processes  such  as  charge-caerier-phonon  coupling  is  crucial 
for  accurate  modeling  of  the  electronic  properties  of  semiconductor  devices.  A  powerful  and 
conceptually  appealing  treatment  of  carrier-phonon  coupling  is  provided  by  the  Feynnan  Path- 
Integral  (FPI)  formalism  of  quantum  mechanics  (1,2|.  The  power  of  the  FPI  formalism  in  treating 
coupling  processes  was  demonstrated  first  in  quantum  electrodynamics  (QED)  (3j,  but  soon  its  ad¬ 
vantages  in  the  formally  similar  treatment  of  carrier-phonon  coupling  in  bulk  semiconductors  also 
was  realised  [4].  Further,  the  reduction  in  device  sizes  to  dimensions  on  the  order  of  the  thermal 
de  Broglie  wavelength  of  charge  carriers  offers  few  conceptual  difficulties  for  the  fully  quantum 
mechanical  FPI  formalism,  while  greatly  reducing  the  value  of  the  semiclassical  treatment  of 
carrier-phonon  scattering  as  spatially  localized,  instantaneous  events.  However,  the  reduction  of 
device  dimensions  does  pose  great  practical  problems  to  application  of  the  FPI  formalism.  For 
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AaSTRACT 

A  hydrodynamic  hot  electron  model  is  used  to  study  electron  transport  through  a  submicron  Af  **“  -  iV  -  iV'*’  GaAf 
structure.  This  study  ts  u«sd  to  investigate  improvements  which  the  unique  features  of  this  model  offer  to  analysis 
2^  devices  operating  under  nonstationary  transport  conditions.  The  model  is  bued  upon  semiclawical  ‘^hydrody* 
aarmc"*  conservation  equations  for  the  average  carrier  density,  momentum  and  energy.  The  general  model  include 
particle  relaxation  times,  momentum  relaxation  times,  energy  relaxation  times,  electron  temperature  tensors  and 
heat  flow  vectors  u  a  function  of  average  carrier  energy  for  the  F,  Y  and  L  valleys  of  GaAs.  For  this  study,  we 
utilized  a  simplified  single  electron  gas  version  of  our  model  to  clearl;.  -eal  the  impact  of  the  nonstationary  terms 
io  the  model.  Results  from  both  a  drift -diffusion  model  approach  a  a  Monte  Carlo  analysis  are  used  to  show 
the  relative  accuracy  and  facility  this  new  model  offers  for  invest;  ing  practical  submicron  device  structures 
operating  under  realistic  conditions. 


KEYWORDS 

Hot-electron;  hydrodynamic  tiansport  model;  Monte  Carlo  method;  nonequilibrium  transport;  nonstationary 
transport;  self-consistent  potentials. 


INTRODUCTION 

Recently,  there  hu  been  a  growing  interest  in  the  use  of  hydrodynamic  conservation  models  to  study  electron 
transport  (Goldsman  and  Frey,  1988;  Satidbom,  Rao  and  Blakey,  1989).  The  impetus  for  choosing  this  approach 
is  very  clear.  Hydrodynamic  models  have  the  ability  to  include  nonstationary  effects  (Hcas  and  lafrate,  1988) 
and  hence  are  superior  in  physical  detail  to  simple  drift-diffusion  models.  Also,  while  Monte  Carlo  methods  can 
euily  incorporate  complicated  band  structure  and  detailed  scattering  rates,  hydrodynamic  models  require  less 
computation  time  to  generate  solutions  and  possess  macroscopic  terms  which  offer  greater  physical  insight  and 
intuition.  However,  there  has  been  some  question  at  to  the  accuracy  of  these  models  u  opposed  to  Monte  Carlo 
methods  due  to  the  inherit  neglect  of  ensemble  effects  (Crandle,  East  and  Blakey,  1989).  In  this  paper,  we  present 
the  results  for  a  transport  study  of  a  submicron  Af  -  AT  -  GaAs  device  structure  to  evaluate  the  importance  of 
the  nonstationary  terms  In  our  new  hydrodynamic  model.  Features  of  this  study  include  :  (1)  The  incorporation 
of  self-consistent  potentials  by  the  simultaneous  solution  of  the  conservation  equations  and  Poisson*!  equation.  (2) 
The  incorporation  of  realistic  boundary  conditions  on  the  device*!  active  region  by  including  an  exact  doping  profile 
for  the  entire  structure  and  applying  boundary  conditions  to  the  contact  points  of  the  Af  ^  regions.  The  multipoint 
nonlinear  boundary  value  problem  theA  results  from  applying  the  transport  model  to  the  device  structure  is  solved 
^sing  an  efficient  local  nonlinear  solver  combined  with  a  perturbation-in-doping  based  continuation  method.  The 
electron  transport  results  are  compared  to  results  from  both  a  drift -diffusion  model  and  a  self-consistent  ensemble 
Monte  Carlo  analysis  of  the  same  GaAs  device  structure. 


1347 


POLAR  OPTICAL  PHONON  SCATTERING  OF 
CHARGE  CARRIERS  IN  ALLOY  SEMICONDUCTORS: 
EFFECTS  OF  PHONON  LOCALIZATION 


L.  F.  Registcr^*\  M.  A.  Littlejohn^®’*^  and  M.  A.  Stroscio^^**^ 


f®^Electricai  and  Computer  Engineering  Department 
North  Carolina  State  University 
Raleigh,  North  Carolina  27695-7911 

Army  Research  Office 

Research  Triangle  Park,  North  Carolina  27709 


ABSTRACT 

the  effects  of  spatial  localization  of  phonons  or  their  correlation  functions  to  finite  distances  in  alloy  semiconduc¬ 
tors  on  polar  optical  phonon  scattering  of  hot  carriers  are  modeled  analytically.  Despite  the  possibility  of  increased 
numbers  of  carrier-scattering-active  phonon  modes,  it  is  demonstrated  that  phonon  localization  in  alloys  should 
have  little  if  any  effect  on  the  total  polar  optical  scattering  rate  for  charge  carriers  coupled  to  equilibrium  phonon 
populations.  Further,  it  is  demonstrated  that  phonon  localization  may  have  a  beneficial  effect  on  hot  carrier 
transport  by  reducing  the  possibility  of  e.xciting  nonequilibrium  phonon  populations.  These  results  are  obtainq^ 
without  assuming  any  specific  functional  form  or  degree  of  phonon  localization;  rather,  calculations  rely  on  the 
inherent  orthogonality  and  mathematical  completeness  of  the  classical  vibrational  modes  over  the  crystal  lattice 
degrees  of  freedom. 

KEYWORDS 

Phonon  localization,  polar  optical  phonon  scattering,  Frohlich  coupi  hot  electron  transport,  alloy  disorder. 


INTRODUCTION 

Alloy  scattering  of  charge  carriers,  the  scattering  of  carriers  due  to  the  aperiodic  arrangement  of  atoms  in  alloy 
semiconductors,  long  hu  been  recognized  as  a  significant  carrier  scattering  process  in  these  materials  (for  example, 
Littlejohn  and  colleagues,  1978  and  references  therein).  More  recently  it  has  been  recognized  that  alloy  disorder' 
similarly  may  affect  phonons:  In  addition  to  either  one-  or  two-mode  behavior  depending  on  alloy  type  and  pos¬ 
sibly  composition  (for  example.  Barker  and  Sievers,  1975),  Raman  light  scattering  studies  of  the  phonon  spectra 
of  alloy  semiconductors  exhibit  broadening  and  asymmetry  in  the  Raman  line  shapes  that  has  been 

escribed  to  spatial  localization  of  the  phonon  modes  or  at  leut  their  correlation  functions  (Jusserand  and  Sapriel, 
1981;  Krabach  and  colleagues,  1983;  Paranthal  and  Poliak,  1984).  In  theory,  as  a  consequence  of  the  Htutnhtr^ 
uncerfainfy  pnncfp/e,  the  localization  of  optical  phonon  states  leads  to  a  breakdown  in  the  usual  q  =  0  Raman 
selection  rule  for  the  phonon  crystal  momentum,  allowing  photon  coupling  to  phonons  with  a  spread  of  nominal 
f  values,  and  with  a  corresponding  observable  spread  in  energy  values.  Parayanthal  and  PoUak  (PP)  (1984), 
assuming  a  Gaussian  form  for  the  localization,  have  calculated  correlation  lengths  based  on  the  degree  of  observed 
broadening  In  the  Raman  lines  of  as  little  as  62  and  80  A  for  the  .4/.4j-Uke  and  GaAs-like  modes  respectively 
in  an  AiotGuc.iAs  alloy  sample,  values  representing  significant  spreads  in  the  crystal  momentum  of  the  Raman 
active  phonon  modes.  These  observable  effects  of  alloy  disorder  in  semiconductors  on  photon-phonon  coupling 
suggest  that  alloy  disorder  also  may  have  a  significant  effect  charge-carrier-phonon  coupling  —  in  addition  to 
representing  a  significant  carrier  scattering  process  unto  itself.  Recently,  however,  Kash,  Jha  and  Tsang  (KJT) 
(1987)  studied  alloy  hot-electron-optical-phonon  coupling  both  experimentally  in  two- mode  alloys 

»nd  one-mode  /n,Gai-,.4j  alloys  using  picosecond  time-domain  pump-probe  Raman  techniques,  and  theoretically 
for  polar  [Frohlich  (Frohlich,  1954))  coupling  in  the  former  allowing  for  variations  in  the  phonon  frequencies  and 
and  polarization  strengths  of  the  the  two  branches  of  the  longitudinal-optical  phonon  spectrum  as  a  function  of 
alloy  composition  x.  In  contrast  to  the  observable  effects  of  alloy  disorder  on  photon-phonon  coupling,  th«f  exper- 
iraental  measurements  of  the  lifetimes  and  sizes  of  small  nonequilibrium  phonon  populations  in  these  alloys  exhibit 
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Intrinsic  and  Extrinsic  Response  of  GaAs  Metal- 
Semiconductor-Metal  Photodetectors 

WACLAW  C.  KOSCIELNIAK.  JEAN-LUC  PELOUARD.  and  MICHAEL  A.  LITTLEJOHN,  senior  member,  ieee 
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Abstract— G»As  mcttl-scmiconductor-meial  (.MSM)  photodeteclors 
have  been  simulated  using  a  self-consistent  Monte  Carlo  (MC)  method. 
Intrinsic  device  properties  are  discussed  in  terms  of  MC  electron  and  hole 
transport  under  low  illumination  intensity.  Parasitic  circuit  elements  are 
then  introduced  to  more  closely  model  realistic  devices  using  the  .MC 
results  in  a  circuit  simulator.  Intrinsic  devices  with  0.5  and  I  spacing 
between  fingers  are  dominated  by  stationary  high-field  transport.  Surpris¬ 
ingly,  full-width-half-maximum  (FWHM)  of  0.5  and  I  fim  detectors  with 
parasitica  is  4.3  and  3.8  ps,  respectively.  However,  the  I  /im  detector 
exhibits  a  long  hole  tail  and  transient  oscillations.  Thus.  FWHM  results 
(and  intrinsic  device  response)  can  be  inadequate  predictors  of  ultimate 
frequency  response  and  scaling  behavior.  However,  an  estimate  of 
maximum  repetition  frequency  gives  =  92  GHz  for  the  0.5  ^m 
device,  consistent  with  experimental  data. 


Recently,  rapid  progress  in  GaAs  metal-semiconduc- 
tor-metal  (MSM)  photodetcctor  performance  has  been 
achieved  [l)-[3].  The  fastest  reported  MSM  detector  has  an 
intrinsic  bandwidth  of  105  GHz  [3],  competitive  with  other 
Schottky  photodiodes  [4],  [5].  However,  the  MSM  photode¬ 
tector  is  better  suited  for  monolithic  integration  [6)-(9]. 

To  design  MSM  photodetectors  with  a  required  bandwidth 
and  responsivity,  one  must  identify  performance-limiting 
factors.  These  factors  include  intrinsic  transport  of  photoexci- 
ted  carriers  between  electrodes  and  transmission  line  effects 
(parasitics)  of  the  contact  fingers.  Ideally,  a  coupled  field 
equation/transport  equation  model  is  needed  to  accurately 
simulate  these  factors.  Such  a  model  is  very  computer 
intensive  and  difficult  to  interpret. 

We  have  used  a  self-consistent  Monte  Carlo  (MC)  model  to 
simulate  intrinsic  transport  and  a  lumped-constant  circuit 
model  for  parasitics.  The  MC  results  are  used  as  an  input  for 
SPICE  simulation  of  the  circuit  model.  This  approach  is 
relatively  easy  to  implement.  However,  the  interpretation  of 
the  results  depends  critically  on  the  parasitic  lumped-constant 
circuit  parameters.  The  purpose  of  this  paper  is  to  present  new 
simulation  results  from  an  MC  study  of  MSM  photodetectors 
with  0.5  and  1  fim  finger  ^Mcings.  A  modified  parasitic  circuit 
gives  good  agreement  with  recent  experimental  data  for  a 
device  with  0.5  /tm  interfinger  spacing. 

The  structure  of  the  photodetector  with  planar  multifinger 
(Al)  contacts  is  shown  in  Fig.  1(a).  Unintentionally  doped 
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(b) 

Fig.  I.  (a)  SchemaiN'  -tructure  ot  a  high-speed  GaAs  MSM  photodetector. 
The  finger  width>  ,nd  iniertinger  spacings  are  either  O.S  or  I  jim.  The 
photodetectors  wu:i  a  O.S  and  I  ftm  interfinger  distance  have  it  and  6 
fingers,  respective .  ib)  Energy  band  diagram  of  photodetectors  at  a  bias 
of  1  V. 

GaAs  (Nd  =  5  <  10‘^  cm*’)  simulates  devices  processed 
either  on  epitaxial  layers  or  semi-insulating  substrates.  The 
contact  barrier  heights  (0.78  eV  and  0.65  eV  for  electrons  and 
holes,  respectively)  are  much  larger  than  kT.  Thus,  current 
injection  from  the  contacts  is  negligible. 

The  MC  algorithm  has  been  previously  described  in  detail 
[10],  A  10  fs  time  step  is  used  in  the  simulations  repotted  here. 
Since  it  is  longer  than  the  electromagnetic  wave  propagation, 
this  value  is  small  enough  for  accurate  numerical  analysis  yet 
long  enough  to  retain  validity  of  Poisson’s  equation.  This  is 
indicated  in  the  calculated  conduction  band  profiles  shown  in 
Fig.  1(b).  The  maximum  density  of  photoexcited  carriers  is 
1.3  X  10'^  cm'  which  prevents  appreciable  fluctuations  of 
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The  linear-chain  approximation  is  used  to  calculate  the  spectrum  of  confined  longitudinal- 
optical  (LO)  phonon  fre  ^uencies  in  short-period  strained-layer  superiattices.  The  frequencies 
for  confined  LO-phonon  modes  are  report^  explicitly  for  the  case  of  GaAs/GaP  short-period 
strained-layer  superiattices  grown  in  the  (001)  direction.  These  results  are  compared  with  the 
few  existing  experimental  measurements  for  such  superiattices. 


I.  INTRODUCTION 

Recently  phonons  in  semiconductor  heterostructures 
have  received  much  attention.  These  efforts  are  due,  in  part, 
to  the  fact  that  scattering  by  the  longitudinal-optical  (LO) 
phonon  modes  is  an  important  energy  loss  mechanism  for  a 
wide  variety  of  III-V  semiconductor  devices.'  Raman  scat¬ 
tering  measurements  have  been  used  to  study  the  properties 
of  phonon  modes  in  quantum  wells  and  superiattices.’'*  The 
binary  GaAs/AlAs  superiattices  have  been  experimentally 
investigated  most  extensively  due  to  the  simplicity  in  analy¬ 
sis  (compared  to  tenary  GaAs/AlGaAs  superiattices)  and 
the  material  familiarity.  In  these  structures,  phonon 
modes  propagating  along  the  direction  normal  to  the  hetero¬ 
junctions  are  observed  to  have  two  different  types  of  charac¬ 
teristics:  The  phonon  modes  in  the  acoustic  branch  are  prop¬ 
agative  across  the  heterojunctions  and,  as  a  result,  are 
“folded”  into  the  reduced  Brillouin  zone,^  while  the  modes 
in  the  optical  branch  are  “localized”  in  each  layer  leading  to 
descriptions  of  confined  modes  or  slab  modes,*  In  addition, 
the  localization  of  phonon  modes  near  the  heterojunctions 
leads  to  the  well-known  interface  modes.*  Furthermore,  it  is 
suggested  that  the  Raman  frequency  shifts  for  confined  LO- 
phonon  modes  can  be  used  as  a  means  to  determine  the  bulk 
phonon  dispersion  relation  (i.e.,  the  equivalent  wave-vector 
model)*  and  as  a  probe  to  characterize  the  quality  of  inter¬ 
faces.*  The  frequency  shift  due  to  confinement  has  been  ana¬ 
lyzed  theoretically  as  well.  The  results  obtained  by  a  simple 
linear-chain  model  show  a  remarkable  agreement  with  ex¬ 
perimental  data  except  in  the  monolayer  superiattices.*  In 
more  detailed  analyses,  the  phonon  dispersion  relation  is  cal¬ 
culated  in  a  plane  parallel  to  the  heterojunction.^'* ' 

The  advent  of  strained-layer  (i.e.,  pseudomorphic) 
structures  adds  an  interesting  effect  on  the  confined  phonon 
modes.  The  mismatch  in  lattice  constants  subjects  these 
structures  to  the  biaxial  stress,  and  results  in  the  strain-in¬ 
duced  shifts  in  optical  phonon  frequencies.  The  measure¬ 
ment  of  Raman  shifts  for  strained-layer  quantum  wells  and 
superiattices  has  been  accomplished  for  a  relatively  small 
number  of  III-V  systems;  most  notably  the  GaSb/AlSb  su¬ 


periattices  by  Jusserand  et  al,'^  and  Santos  et  al,,'^  the 
GaAs/InAJAs  superiattices  by  Nakayama  et  al,'*  and  the 
GaAs/GaP  superiattices  by  Armelles  et  a/.'*  These  experi¬ 
ments  clearly  exhibit  the  Raman  shifts  due  to  both  strain  and 
confinement.  In  the  analysis  of  experimental  data,  the  con¬ 
tributions  by  these  two  effects  are  in  general  treated  sepa¬ 
rately  considering  one  effect  at  a  time.  This  approach  is  read¬ 
ily  justifiable  when  the  mismatch  is  not  large  sich  as  in 
GaSb/ AlSb  systems  ( less  than  1  % ) .  However,  as  the  strain 
becomes  sig^' ‘leant  (such  as  in  GaAs/InAs  systems  where 
the  mismatL  >$  approximately  7%),  the  coupling  between 
two  effects  ly  becomes  non-negligible.  It  is,  thus,  neces¬ 
sary  to  stud;  he  importance  of  the  coupling  effect,  and  de¬ 
velop  a  modi  which  can  directly  predict  phonon  frequencies 
in  this  situation. 

In  this  paper  we  present  a  model  for  determining  the 
frequencies  of  confined  LO  phonons  in  strained-layer  III-V 
superiattices  grown  in  the  ((301  >  direction.  A  discrete  (rath¬ 
er  than  continuum )  approach  based  on  a  linear-chain  model 
is  applied.  The  frequency  shift  arising  simultaneously  from 
both  confinement  and  strain  is  accounted  for  by  properly 
modifying  the  force  constants  to  incorporate  the  strain  ef¬ 
fects  into  the  linear-chain  model.  For  a  specific  example,  the 
frequency  shifts  are  calculated  for  GaAs/GaP  short-period 
superiattices,  and  compared  with  the  existing  Raman  data  of 
Armelles  et  a/.'*  for  the  lowest  order  confined  LO  modes. 
Our  model  will  be  of  importance  in  the  strained-layer  struc¬ 
tures  when  the  strains  have  to  be  estimated  rather  accurate¬ 
ly:  The  diagnosis  of  heterojunction  interface  quality*  is  one 
example,  and  the  study  of  carrier  transport  is  another.  Since 
the  strains  in  individual  layers  modify  the  intrinsic  physical 
properties  such  as  band  structures,'*  the  strain-induced  ef¬ 
fects  can  be  applied  to  tailor  the  device  structures  for  a  possi¬ 
ble  reduction  of  carrier- LO-phonon  interactions.'’ 

II.  MODEL 

Figure  1  shows  a  typical  schematic  drawing  for  the  lin¬ 
ear-chain  model  in  a  superlattice.  As  can  be  seen  in  the  fig¬ 
ure,  the  z  axis  is  chosen  as  the  growth  direction  of  the  lattice 
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Abstract— Wt  present  calculations  of  the  electron  und  hole  compo¬ 
nents  of  dark  current  in  a  GaAs  metal-semiconductor-metal  iMS.M) 
photodetector.  .A  quantum-mechanical  model  is  deseloped  that  de¬ 
scribes  the  electron  and  hole  transport  behavior  in  the  contact  reuions 
which  is  used  to  determine  dark  current  as  a  function  of  electric  Held. 
The  model  reduces  to  a  conventional  thermionic  emission  model  if  an 
ideal  barrier  transmission  coefficient  is  assumed.  In  order  to  assess  the 
accuracy  of  the  model,  photodetectors  have  been  fabricated  and  tested. 
Theoretical  calculations  and  experimental  data  are  compared  and  pood 
agreement  is  obtained.  Possible  modifications  to  enhance  the  usefulness 
of  the  model  are  discussed. 


I.  I.NTRODUCTION 

ALLIUM  ARSENIDE  (GaAs)  metal-semiconduc- 
tor-metal  (MSM)  photodetectors  have  recently  been 
used  to  detect  ultrashort  optical  pulses  [If  With  a  planar, 
interdigitated  pattern  delineated  by  electron-beam  lithog¬ 
raphy.  an  intrinsic  detector  bandwidth  in  excess  of  100 
GHz  has  been  achieved  {2).  Increased  research  in  this  area 
is  stimulated  by  possible  development  of  ultra-high-speed 
monolithic  optoelectronic  integrated  circuits  (lOEC's).  in 
which  the  MSM  photodetector  would  serve  as  an  input 
device  for  optical  signals  and  would  be  followed  by  a  sig¬ 
nal  processing  circuit.  Applications  of  such  circuits  for 
telecommunications  and  optical  data  processing  are  fea¬ 
sible  in  the  near  future,  especially  since  monolithic  inte¬ 
grated  amplifiers  operating  up  to  60  GHz  have  already 
been  demonstrated  [3]. 

To  better  understand  major  factors  limiting  speed  per¬ 
formance  of  this  detector,  detailed  numerical  simulations 
have  been  undertaken  by  the  present  authors  (4|-(61  and 
others  (2).  In  order  to  achieve  a  large  signal-to-noise  ra¬ 
tio.  detectors  with  low  dark  .current  and  low  noise  are  re¬ 
quired.  However,  only  basic  aspects  of  dark  current  of  the 
MSM  photodetectors  have  been  addressed  (7]-(91.  In  the 
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following,  we  discuss  a  one-dimensional  model  for  deter¬ 
mining  the  electron  and  hole  components  of  the  dark  cur¬ 
rent.  This  approach  integrates  solutions  of  the  time-in¬ 
dependent  Schrodinger  equation  with  thermionic  emission 
models.  The  calculations  are  compared  with  data  from 
fabricated  devices,  and  good  agreement  between  theory 
and  experiment  is  obtained.  Also,  suggestions  for  further 
extension  of  the  model  are  given. 

11.  Device  Fabrication 

A  schematic  diagram  of  the  experimental  photodetector 
is  shown  in  Fig.  1 .  The  multifinger  pattern  was  fabricated 
on  a  high-purity,  undoped  substrate.  A  standard  lift-off 
process  (chlorobenzene-soaked  positive  photoresist)  was 
used  to  delineate  4000-A-thick  A1  fingers.  Shortly*before 
A1  deposition,  the  native  oxides  were  removed  by  immer¬ 
sion  of  the  sai’  les  in  1  HCl :  1  H.O.  A  recess,  approxi¬ 
mately  600  A  op.  was  etched  to  remove  any  photoresist 
residues  betwe.  i  fingers.  Etching  was  performed  in  a  so¬ 
lution  of  3  H:l  O4;  1  H.OiiSO  H;0  followed  by  a  short 
immersion  in  I  .  (Cl:  1  HyO.  Our  experiments  indicate  that 
the  recess  is  likely  to  be  responsible  for  a  relatively  high 
breakdown  voltage  of  80-120  V  obtained  in  the  best  de¬ 
vices.  Experimental  results  are  presented  for  photodetec¬ 
tors  with  3-5-^4m  distance  between  fingers  and  an  active 
area  of  83  /im  x  83  /im  or  100  /xm  x  100  /xm.  suitable 
for  coupling  to  multimode  optical  fibers. 

Simultaneously  with  the  photodetectors.  A1  Schottky 
diodes  (squares.  7.5  x  10"'  cm’)  were  fabricated  on 
Si-doped  (2  x  lO’’  cm"')  substrates.  Earlier.  Au- 
GeNi/Au  ohmic  contacts  were  alloyed  to  the  back  sides 
of  the  n-type  samples.  The  barrier  height  (0„  =  0.75- 
0.77  eV)  and  the  ideality  factor  {n  =  1.04-1.05)  were 
determined  from  l-V  characteristics  of  these  diodes.  The 
near-unity  ideality  factor  indicates  a  high-quality  Al- 
GaAs  interface  and  absence  of  native  interfacial  oxides. 
This  information  is  important  since  we  neglect  the  native 
oxides  in  the  model  development.  In  addition,  the  change 
in  /-P' characteristics  upon  various  annealing  temperature 
cycles  up  to  470®C  was  measured  in  order  to  obtain  ad¬ 
ditional  information  about  the  nature  of  the  metal-semi- 
conductor  interface.  There  is  an  initial  increase  in  barrier 
height  for  temperatures  up  to  250®C,  followed  by  a  slow 
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A  linear-chain  model  wuh  variable  force  constants  for  the  couplings  at  heicrojunc- 
tions  IS  used  to  calculate  the  spectrum  of  confined  loncinidinal-optical  (LO)  phonon  fre¬ 
quencies  in  a  shon-penod  strained-layer  supcrlanicc.  Even  for  the  case  of  a  supcrlarace 
with  only  rwo  atomic  monolayers  in  each  superlaracc  layer,  it  is  dcmonstraicd  that  the 
frequencies  of  contined  LO-phonon  modes  are  only  weakly  dependent  on  variations  in 
interfacial  force  constants. 


1.  Introduction 

Phonons  In  semiconductor  hcterostructurcs  have 
received  much  attention.  These  efforts  are  due.  in  pan.  to 
the  fact  that  scattering  by  the  longitudinal-optical  (LO) 
phonon  modes  is  an  important  energy  loss  mechanism  for 
electrons  in  a  wide  variety  of  IH-V  semiconductor  dev¬ 
ices.^  In  addition,  Raman  scattering  measurements  have 
been  applied  to  study  the  properties  of  phonon  modes  in 
quantum  wells  and  superlatticcs.*"^  Recently.  Raman 
spectra  have  been  used  to  characterize  monolayer  varia- 
tiona  in  the  thicknesses  of  quantum  wells  ®  In  the  analysis 
and  interpretation  of  Raman  data  for  confined  LO-phonon 
modes,  it  is  convenient  to  model  the  confined  LO  modes 
as  modes  in  a  uniform  isolated  dielectric  siab.^~‘^  Such 
i  an  approximate  model  Is  based  on  the  as.sumption  that 

j  atomic  force  constants  at  hcterojunction  interfaces  are 

}  identical  with  those  of  the  bulk  or  of  uniform  pseu- 

;  domorphic  layers.  In  this  paper,  the  effect  of  varying 

]  force  constants  at  the  hcterojunction  interfaces  of  a 

J  stramcd-layer,  short-period  superlattice  is  analyzed  by 

^  using  a  linear-chain  model  to  determine  the  corresponding 

:  changes  in  the  frequencies  of  the  confined  LO-phonon 

.  modes.  In  the  studies  reported  hereui.  a  superlattice  with 

only  two  atomic  monolayers  m  each  superlattice  layer  is 
I  modeled  in  order  to  dctemtuie  the  unportance  of  perturba¬ 

tions  in  heterojunction  interface  force  constants  for  a  sys¬ 


tem  where  such  effects  should  be  large.  In  particular,  fre¬ 
quencies  of  confined  LO-phonon  modes  are  determined 
for  a  GaAs/GaP  short-period  superlanice  with  two  mono- 
layers  per  superlattice  layer.  It  is  found  that  varying  inter- 
facial  torce  constants  by  values  as  extreme  as  10%  results 
m  only  about  a  2%  change  in  the  frequencies  of  the 
confined  LO-phonon  modes.  As  a  practical  matter, 
changes  m  the  frequencies  of  confined  LO-phonon  modes 
will  generally  be  considerably  less  than  2%  since  in  most 
superLittices  and  quantum  wells  the  ratio  of  the  number  of 
bonds  ai  inicrfaccs  to  the  number  of  bonds  removed  from 
interfaces  is  less  than  for  the  two-monolayer  per  laver  case 
exammed  herein. 

2.  Linear-Chain  Model 
with  Variable  Interfacial  Force  Coiistanis 

to  order  to  introduce  the  basic  properties  of  the 
lincar-cham  mode!  used  in  this  analysis,  we  will  consider 
mitiailv  the  case  where  interface  effects  are  absent. 
Figure  1  depicts  a  schematic  for  the  itncar-cham  system 
used  (0  study  ihc  role  of  force-constant  variations  at 
hctcroiunciion  mtcrfaccs.  The  z  axis  is  chosen  as  the 
growth  direction  ot  the  lattice.  Hence,  for  the  longitudinal 
modes  ot  interest,  the  atomic  displacements  arc  in  the  z 
direction  In  the  absence  of  strain,  the  force  constant  m 
each  laver  is  estunated  based  on  the  frequencies  of  the 


4l 


0749-6036/90/0201 15-04  S02.00;0 


1 990  Academic  Press  Limited 


Appendix  C 

Resumes  for  Principal  Investigators,  Dr.  M.  A.  Littlejohn  and  Dr.  K.  W.  Kim. 


BIOGRAPHICAL  SUMMARY  FOR  MICHAEL  A.  LITTLEJOHN 


BiisingM  Addresa 

Department  ctf  Electrical  and  Computer  Engineering 
North  Carolina  State  Univenity 
Box  7911 

Raleigh,  North  CaroUna  27695-7911 
(919)737-2336 


Education 

North  Carolina  State  University 
North  (Carolina  State  University 
North  Carolina  State  University 


1962 

1964 

1967 


B.S.  Electrical  En^neering 

M.S.  Electrical  Engineering 

Ph.D.  Electrical  En^eering 


Industrial  aniL 
1967-prescnt 


1977-1981, 

198S-present 

1981-1983 

1984-1985 

1988-1989 


Academic  Experience 

Assistant  I^fessor/Associate  Professor/Professor,  Department  of  Electrical 
and  Gxnputer  Engineering,  Nordi  (^lina  State  University,  Raleigh, 
North  (Molina. 

Technical  Staff  Member,  U.S.  Army  Research  Office,  Research 
Triangle  Park,  North  (Carolina. 

Director,  Microelectronics  Programs,  School  of  Engineering,  North 
Carolina  State  University,  Raleigh,  North  Carolina. 

Associate  Dean  of  En^eering  fw  Research  Ptograms,  North  (Carolina 
State  University,  Raleigh,  North  (Carolina. 

Associate  Director,  NSF  Engineering  Research  Onter  On  Advanced 
Electronic  Materials  Processing. 


Prgfcssiffnal  Attujtffti 

Member  at  Large,  Executive  (Committee  of  the  Electronics  Division  the  Electrochemical 

Society,  1983-1988. 

Member,  National  Academy  of  Sciences  National  Research  (Council  Committee  on 
Recommendations  for  U.S.  Army  Basic  Scientific  Research,  1985-1986. 

Member,  Technical  Advisory  Board,  IEEE  En^neering  Research  and  Development 
CkMomittee,  1986*1989. 

Associate  Editor,  Electronics  Division,  Journal  of  the  Electrochemical  Society,  1987-1990. 

Member,  NSF  Committee  on  Research  Trends  and  Opi)ortunities  for  the  Metallurgy, 
Polymers,  and  Ceramics  Section,  Division  of  Materials  Research,  1983-1985. 

Organizer  and  Program  (Chairman,  NSF  Wtxkshop  on  The  Future  of  Microstructure 
Technology,  1985. 

Cb-organizer,  NSF  Workshop  on  Co(^>eration  and  Sharing  Among  U.S.  Microelectronics 
Centers,  1983. 


Society  Memberships  and  Honors 

Institute  of  Electrical  and  Electronics  Engineers 
American  Soci^  for  Engineering  Education 
American  PhyriM  Society 
Electrochemi^  Soc^ 

Certification  of  Achievement  Fix'  Patriotic  Civilian  Service,  U.S.  Army,  1989 

RJ.  Reynolds  Industries,  Inc.,  Award  for  Excellence  in  Teaching  and  Research,  1983-1987 

Alcoa  Foundation  Distinguished  Research  Award,  1983 

N.C.  State  University  Alumni  Award — University  Distinguished  Alumni  Professor,  1980 
Western  Electric — ASEE  Fund  Award  for  Excellence  in  Teaching  and  Research  (1978) 
Sigma  Xi  Outstanding  Young  Scientist  Award,  1976 
Eia  Kappa  Nu  Outstanding  Teacher  Award,  1974  &  1975 


Field  of  Research  Interest 

Seimconductor  device  simulation  and  modeling,  1II>V  compound  semiconductor  materials  and 
devices,  hot  electron  tran^tt  in  semiconductors,  ion  implantadon  and  radiation  damage  in 
semiconductors,  thin  films  and  oxides  on  semiconductors,  defects  in  semiconductors. 

Ten  Most  Recent  Refereed  Publications  (From  a  total  of  931 

1 .  W.C  Koscielniak,  J.L.  Pelouard  and  M.A.  Littlejohn,  “Intrinsic  and  Extrinsic  Response  of 
GaAs  Metal-Semiconductor-Metal  Photodetectors,”  Photonics  Tech.  Len.,  2, 125, 1990. 

2.  LF.  Register,  M.A.  Littlejdm  and  MA.  Stioscio,  “Polar  Optical  Phonon  Scattering  of 
Charge  Carriers  in  Alloy  ^miconductors:  Effects  of  Phonon  Localization,”  Solid-State 
Election.,  22. 1387,  1989. 

3.  D.L.  Woolard,  JX.  Pelouard,  RJ.  Tiew,  MA.  Littlejohn  and  C.T.  Kelley,  “Hydrodynamic 
Hot-Electron  Transport  Simulation  Based  on  the  Monte  Cario  Method,”  Solid-State  Electron., 
22.  1347.  1989. 

4.  W.C.  Koscielniak,  M.A.  Littlejohn  and  J.L  Pelouard,  “Analysis  of  a  High-Speed  GaAs 
Metal-Semiconductor-Metal  (MSM)  Photodetector  with  0.1  pm  Finger  Spacing,”  THEE 
Election  Dev.  Lett.,  10, 221  (1989). 

5.  W.C  Koscielniak,  J.L.  Pelouard  and  M.A.  Littlejohi^  “Dynamic  Behavitv  of  Photocuirent  in 
a  GaAs  Metal-Sen^conductor-Metal  Photodetector  with  Sub-Half  Micron  Electrode  Pattern,” 
Appl.  Phys.  Lett,  2X  567  (1989). 

6.  K  J.  Brennan,  D.R  Park,  K.  Hess  and  M.A.  Littlejdin,  “Theory  of  the  Velocity  Field 
Relation  in  AlGaAs,”  Jour.  Appl.  Phys.,  SL  5004  (1988). 

7.  D.L.  Woolard,  RJ.  Trew  and  M.A.  Littlejohn,  “Hydrodynamic  Hot  Election  TranspcM 
Model  With  Monte  Carlo-Generated  Transport  Parameters,”  Solid  State  Electronics,  3).  571 
(1988). 

8.  LF.  Register,  M.A.  Littlejohn  and  MA  Stioscio,  “Feynman  Path  Integral  Study  of  (Confined 

C^ers  Subject  to  a  Statistical  Potential,”  Solid  State  Electronics,  2L  563  (1988).  * 

9.  Li^.  Register,  MA  Stioscio  and  M.A.  Littlejohn,  “Numerical  Evaluation  of  the  Feynman 
Integral-(Dver-Paths  in  Real  and  Imaginary  Time,”  Supeilattices  and  Microstiuctures,  i,  61 
(1988). 

10.  C.K.  ^liams,  M.A.  Littlejohn,  T.H.  Glisson  and  J.R.  Hauser,  “Monte  Carlo  Simulation  of 
the  Hall  Effect  in  Degenerate  GaAs,”  Superlattices  and  Microstructures.  22.  P-  725, 1986. 

Associations  With  Graduate  Students/Post-Doctoral  Scholars  For  Last  Five  Years 

(From  a  total  of  57  graduate  students  and  7  post«doctoral$  sponsored) 

1.  H.P.  Belgal,  M.S.E.E..  1988. 

2.  L.F.  register,  Ph.D.,  1989. 

3.  D.L.  Woolaid,  Ph.D,  1990. 

4.  W.C.  Koscielniak,  Ph.D,  1990. 

5.  K.J.  Kim,  Ph.D. 

6.  J.E.  Mona^ian,  I%.D. 

7.  J.L  Pelouard,  Visiting  Assistant  Professor,  1988-1989. 

8.  H.H.  Tian,  Visiting  Assistant  Professor,  1989-1990. 

Scientific  Collaborators  During  Last  Four  Years 

1 .  Dr.  M.A.  Stroscio,  U.S.  Army  Research  (Office,  RTP,  NC 

2.  Dr.  KP.  Brennan,  Georgia  Institute  of  Technology,  Atlanta,  GA 

3.  Dr.  K.  Hess,  University  of  Illinois,  Urbana,  IL 


Ki  Wook  Kim 


Dept,  of  Electrical  and  Con^uter  Engineering 
Box  7911 

North  Carolina  State  University 
Raleigh,  NC  27695-7911 

EDUCATION 

08/88  Ph.D.  in  Electrical  Engineering 

University  of  Dlmois,  Urbana,  Illinois 
Thesis  Title:  Monte  Carlo  Studies  of  Nonlinear  Electron  Transport 
in  ni-V  Semiconductors 

05/85  M.S.  in  Electrical  Engineering 

University  of  Illinois,  Urbana,  Illinois 

02/83  B.S.  in  Electronics  Engineering 

Seoul  National  University,  Seoul,  Korea 

WORK  EXPERIENCE 


Phone:  (919)737-5229 
Birthdate:  September  14, 1960 


08/88  -  present:  Assistant  Professor  (Visiting),  Department  of  Electrical  and  Computer 
Engineering,  North  Carolina  State  University,  Raleigh,  North  Carolina 

06/85  -  08/88:  Graduate  Research  Assistant,  Coordinated  Science  Laboratory,  University  of 

Illinois,  Urbana,  Illinois 


06/83  -  05/85:  Graduate  Research  Assistant,  Department  of  Electrical  Engineering,  Univer¬ 

sity  of  Illinois,  Urbana,  Illinois 


ACADEMIC  HONORS 


IBM  Predoctoral  Fellowship  (1986-1988) 
RESEARCH  INTERESTS 


Semiconductor  physics  and  modeling  of  electronic  and  optoelectronic  devices,  carrier  transport 
in  bulk  and  heterostructures,  quantum  effects,  Monte  Carlo  simulation. 


PUBLICATIONS 


Refereed  Articles  in  Journals 

1.  K.  Kim,  K.  Kahen,  J.P.  Leburton,  and  K.  Hess,  "Band-structure  dependence  of  impact  ion¬ 
ization  rate  in  GaAs,"  J.  Appl.  Phys.  59, 2595  (1986). 

2.  K.  Kim  and  K.  Hess,  "Simulations  of  electron  impact  ionization  rate  in  GaAs  m  nonuni¬ 
form  electric  fields,"  J.  Appl.  Phys.  60, 2626  (1986). 

3.  D.  Arnold,  K.  Kim,  and  K.  Hess,  "Effects  of  field  fluctuation  on  impact  ionization  rates  in 
semiconductor  devices  due  to  the  discreteness  and  distribution  of  dopants,"  J.  Appl.  Phys. 
61,1456(1987). 

4.  K.  Kim,  K.  Hess,  and  F.  Capasso,  "New  effects  of  stmcture  in  momentum  and  real  space 
on  nonlinear  transport  across  heterojunction  band  discontinuities,"  Appl.  Phys.  Lett.  51, 
508  (1987). 


5.  K.  Kim,  B.  Mason,  and  K.  Hess,  "Inclusion  of  collision  broadening  in  semiconductor  elec-» 
tron  transport  simdations,"  Phys.  Rev.  B36, 6547  (1987). 

6.  K.  Kim  and  K.  Hess,  "Ensemble  Monte  Carlo  simulations  of  semiclassical  nonlinear  elec¬ 
tron  transport  across  heterojunction  band  discontinuities,"  Solid-State  Electron.  31,  877 
(1988). 

7.  K.  Kim,  K.  Hess,  and  F.  Capasso,  "A  Monte  Carlo  study  of  electron  heating  and  enhanced 
thermionic  emission  by  hot  phonons  in  heterolayers,"  Appl.  Phys.  Lett.  52, 1167  (1988). 

8.  M.  Kuzuhara,  K.  Kim,  D.  Arnold,  and  K.  Hess,  "Ballistic  electron  transport  across  collec¬ 
tor  barriers  in  AlGaAs/GaAs  hot-electron  transistors,"  Appl.  Phys.  Lett.  52, 1252  (1988). 

9.  K.  Kim  and  K.  Hess,  "Electron  transport  in  AlGaAs/GaAs  tunneling  hot  electron  transfer 
amplifiers,"  J.  Appl.  Phys.  64, 3057  (1988). 

10.  J.  M.  Higman,  K.  Kim,  K.  Hess,  T.  van  Zutphen,  and  H.  M.  J.  Boots,  "Monte  Carlo  simu¬ 
lation  of  Si  and  GaAs  avalanche  electron  emitting  diodes,"  J.  Appl.  Phys.  65,  1384 
(1989). 

11.  M.  Kuzuhara,  K.  Kim,  and  K.  Hess,  "Transient  simulation  of  AlGaAs/GaAs/AIGaAs  and 
AlGaAs/InGaAs/AlGaAs  hot-electron  transistors,"  IEEE  Trans.  Electron  Devices  ED-36, 
118(1989). 


12.  K.  W.  Kim,  M.  A.  Stroscio,  and  J.  C.  Hall,  "Frequencies  of  confined  longitudinal-optical 
phonon  modes  in  GaAs/GaP  short-period  strained-layer  superlattices,"  J.  Appl.  Phys.  67, 
6179  (1990). 

13.  M.  A.  Stroscio,  K.  W.  Kim,  and  J.  C.  Hall,  "Variation  in  frequencies  of  confined 
longitudinal-optical  phonon  modes  due  to  changes  in  the  effective  force  constants  near 
heterojunction  interfaces,"  Superlatt.  Microstruct.  7, 115  (1990). 

14.  M.  A.  Stroscio,  K.  W.  Kim,  M.  A.  Littlejohn,  and  H.  Chuang,  "Polarization  Eigenvectors 
of  Surface-Optical-Phonon  Modes  in  a  Rectangular  Quantum  Wire,"  Phys.  Rev.  B.  42, 
1488  (1990). 

15.  K.  W.  Kim  and  M.  A.  Stroscio,  "Electron-Optical-Phonon  Interaction  in  Binary/Temary 
Heterostructures,"  submitted  to  Appl  Phys.  Lett. 

16.  M.  A.  Stroscio,  K.  W.  Kim,  and  S.  Rudin,  "Scattering  Rate  for  the  Electron-LO-Phonon 
Interaction  in  Polar  Semiconductor  Quantiun  Wires,"  accepted  for  publication  in  Phys. 
Rev.  B. 

17.  K.  W.  Kim,  H.  Tian,  and  M.  A.  Littlejohn,  "Analysis  of  Delta-Doped  and  Uniformly* 
Doped  HEMTs  by  Ensemble  Monte  Carlo  Simulation,"  submitted  to  IEEE  Electron  Dev¬ 
ice  Lett. 

18.  K.  W.  Kim,  M.  A.  Stroscio,  and  A.  Bhatt,  "Electron-Optical-Phonon  Scattering  Rates  in  a 
Rectangular  Semiconductor  Quantum  Wke,"  to  be  submitted  to  Phys.  Rev.  B. 

19.  H.  Tian,  K.  W.  Kim,  and  M.  A.  Littlejohn,  "Influence  of  DX  Centers  and  Surface  States  on 
Delta-Doped  HEMT  Performance,"  to  be  submitted  to  J.  Appl  Phys. 

Refereed  Articles  in  Proceedings  of  Conferences 

1.  K.  Kim,  K.  Hess,  and  F.  Capasso,  "Simulations  of  nonlinear  transport  in  AlGaAs/GaAs 
single  well  heterostructures,"  presented  at  the  5th  Int.  Conf.  on  Hot  Carriers  in  Semicon¬ 
ductors  (July,  1987,  Boston,  Mass.),  published  in  Solid-State  Electron.  31, 349  (1988). 

2.  K.  W.  Kim,  M.  A.  Stroscio,  and  J.  C.  Hall,  "Frequencies  of  Confined  Longitudinal-Optical 
Phonon  Modes  in  Short-Period  Strained  Semiconductor  Superlattices,"  presented  at  the 
SPIE  Int.  Symp.  on  Optical  &  Optoelectronic  Applied  Science  and  Engineering  (July, 
1990,  San  Diego,  Calif.),  Proceedings  to  be  published. 

3.  M.  A.  Stroscio,  K.  W.  Kim,  and  M.  A.  Littlejohn,  "Theory  of  Optical-Phonon  Interactions 
in  Quantum  Wells  and  Quantum  Wires  (Invited),"  to  be  presented  at  the  SPIE  Int.  Conf. 
on  Physical  Concepts  of  Materials  for  Novel  Optoelectronic  Device  Applications 
(October,  1990,  Aachen,  Germany),  Proceedings  to  be  published. 


